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EXECUTIVE SUMMARY 


The design of a modern military force structure is a complicated process. This 
process involves many competing elements. Political pressures, budget constraints and 
foreign threats are some of the main effects in force structuring design. Currently, there 
are several different methodologies for determining the force structure of a naval unit. 
This thesis proposes an alternative methodology for determining the force structure of the 
Turkish Naval Surface Group (TNSG). 

The TNSG force structuring problem (selecting constituent ship types for this 
surface group) can be modeled as a project selection problem. Project selection models are 
examined, and it is considered that the model must be capable of solving multi-objective 
problems (decreasing the cost and increasing the effectiveness of the force for our study). 

Two main models are formulated to solve this force structuring problem. Model 1 
is formulated to determine the number and mix of ship types to purchase or build given a 
constant budget and weighted preferences for each ship type. In this model, the main goal 
is to maximize the effectiveness of the force mtx when a constant budget is given. Model 
2’s formulation uses the weapon requirements needed to meet the threat’s specific 
capabilities (defined in the mission areas of air, surface, submarine, and amphibious 
operations), ship weapon capacities for these requirements, and the weighted preferences 
for the ship types. In second model, the main goal is to minimize the cost of providing a 
fixed force effectiveness. 

There are additional factors, aside from cost and weapon requirements, which 
must be included when considering the mix and size of the naval surface group. These 
factors should include not only the ship attributes such as speed, cost, warfare capabilities, 
and fuel consumption, but with respect to these attributes, the relative advantage held by 
each ship in accomplishing the mission. A survey is a method to obtain subjective 


judgments rating ships with respect to the above factors. The survey judgments are 


XVII 


converted into numeric values, using the Analytic Hierarchy Process which are then used 
as objective function coefficient weights in our model formulations. 

The scenario on which weapon requirements are based 1s a conventional war 
situation with a specific country. Exact values for future weapon requirements will be 
determined by the Turkish Naval Surface Group Staff after examining the force structure 
and weapon capabilities of the threat countries. In this study, however, only approximate 
values are used for requirements and budgets. Hence, this model is a generic model and 
can be used by any country by changing these requirement and budget values. 

This study also provides a sensitivity analysis of the changes in the force mix with 
changes in the future weapon requirements and budget values. It also provides a 
comparison of the models that are used and selects the best model for determining the best 
mix of ship types. 

Several computer packages and programs are used to solve these force structuring 
problems. With the use of software and computer support, it 1s easy to make changes in 
the model values. The optimum ship mixes for the models are found using these computer 
programs. The comparison of the models is made according to the smallest goal deviation 
and the highest benefit-cost ratio, and it is found that Model 2 (minimum cost, fixed 


weapon requirements) gives the best ship mix for TNSG. 
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I. INTRODUCTION 


A. BACKGROUND 


The design of a modern military force structure is a complicated process. This 
process involves many competing elements. Currently, force structure design is affected by 
political pressures, budget constraints and foreign threats. This process impacts on 
Turkey’s ability to protect its national interests in peacetime and war. The main purpose of 
this study is to propose a model which can generate a solution to the Turkish Naval 
Surface Group’s (TNSG) force structuring problem. 

The mission of TNSG is to conduct Anti Surface Warfare (ASUW) operations, 
Anti Submarine Warfare (ASW) operations, Anti Air Warfare (AAW) operations, and 
Amphibious Warfare (AMPW) operations in peacetime (for training) and in war. The 
current TNSG inventory contains only Guided Missile Destroyers (DDG), Guided Missile 
Frigates (FFG), and Guided Missile Fast Patrol Boats (PBFA). However there have been 
some proposals to add other types of surface ships (Cruisers and Corvettes) to the 
inventory. 

Modern cruisers (CG) and guided missile destroyers (DDG) operate in support of 
carrier battle groups, surface action groups, amphibious groups, and replenishment 
groups. Both CGs and DDGs are multi-mission (ASW, AAW, and ASUW) surface 
combatants; however, the CG is also capable of operating independently and is capable of 
serving as a flagship of a surface action group. While guided missile frigates (FFG) can 
fulfill the Protection of Shipping (POS) mission as ASW combatants for amphibious 
expeditionary forces, underway replenishment groups and merchant convoys, they can 
also perform AAW and ASUW missions. But their weapon capacity is limited when 
compared with CGs and DDGs. Figure 1.1 shows an Oliver Hazard Perry class frigate that 
the Turkish Navy is planning to buy from U.S. Navy. Corvettes (FSG) are also multi- 
mission (ASW and ASUW) surface combatants. New types of corvettes also have AAW 


capabilities, with the addition of surface to air guided missiles (SAM). Guided missile fast 


patrol boats (PBFA) are single mission (ASUW) surface combatants. Because of their fast 
speed and addition of the surface to surface guided missiles (SAM), PBFAs are more 
effective in surface warfare. PBFAs also operate in coastal patrol and interdiction 


surveillance. 
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Figure 1.1 Guided Missile Frigate (FFG) . 


Important factors in the TNSG’s ability to conduct their own missions and to assist 
other forces in conducting operations are the ship-type characteristics of speed, ASUW 
capabilities (both open sea and shoreline), ASW capabilities; AAW and Close-in Weapon 
Systems (CIWS) capabilities (against guided missiles and aircraft), AMPW capabilities, 
and surveillance (target reporting) capabilities. These characteristics change according to 
the ship types utilized for a mission. Table 1.1 summarizes surface ship characteristics 
which will be used in our proposed model. Ship characteristics are taken from Jane’s 
Fighting Ships 1996/1997 edition [Ref. 1]. Ship cost is certainly an important 
consideration to take into account when planning within a specified budget. The cost 


values shown in Table 1.1 are taken from Healy [Ref. 2]. 
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Table 1.1 Surface Ships Characteristics. 


B. OBJECTIVE 


The TNSG inventory contains old ships that have been bought from other countries 
and new modern ships. The main goal is to modernize all of the ships in TNSG. Current 
modernization progress 1s mainly focused on building FFGs and PBFAs. The objective of 
this thesis is to provide an alternative methodology for determining the force structure of 
the Turkish Naval Surface Group in this modernization effort. The problem is to determine 
the best mix of ship types given some constraints, such as budgets and/or specific weapon 


requirements. 


The model proposed here obtains survey judgments which rate each surface ship 
type’s relative superiority over the others with respect to selected characteristics (e.g., 
speed, warfare capabilities, surveillance capabilities, and fuel consumption) which enable 
the ship’s operations in a warfare area. The Analytic Hierarchy Process (AHP) is used to 
convert these subjective judgments into numeric values that are used as objective function 
coefficient weights in several goal programming model formulations. The object of each 


formulation is to select the best mix of surface ships by minimizing the total deviation from 


one or more force level goals given certain system constraints, such as budgetary 


requirements, weapon requirements, and/or existing force levels. 
Cc: SCOPE 


The scenario on which weapon requirements are based is a conventional war 
situation with a specific country. Exact values for future weapon requirements will be 
determined by the Turkish Naval Surface Group Staff by examining the force structure 
and weapon capabilities of the threat countries. In this study, however, only approximate 
values are used for requirements and budgets. Hence, this model is a generic model and 
can be used by any country by changing these requirements and budgets. 

This paper 1s organized to logically discuss the modeling process and results. 
Chapter II will consist of a search for the most appropriate project selection model among 
all such models currently used in the operations research literature. Chapter III will 
contain the model development process, and Chapter IV will examine the formulation of 
the models being implemented. Chapter V will present the computer programs used to 
exercise the model as well as model results and output. Chapter VI will contain the 
sensitivity analysis of the models in various situations. The final chapter will present 


conclusions and recommendations. 


Hi. LITERATURE REVIEW 


The force structuring of the Turkish Naval Surface Group (selecting constituent 
ship types) can be viewed as a project selection problem. The intent of this chapter is to 
provide a brief overview and assessment of the literature in the area of project selection. 
The four basic types of project selection models that will be discussed below are subjective 
models, financial analysis models, risk analysis models, and mathematical programming 


models. 
A. SUBJECTIVE MODELS 


The subjective models are the simplest form of formal R&D project evaluation. 
The checklist and scoring models are the subjective models that are used most frequently. 
Liberatore and Titus [Ref. 3] reported the results of an empirical study on the usage of 
quantitative techniques for R&D project management in 40 respondents from 29 random 
Fortune 500 industrial firms. They found that almost half of the respondents had used 


checklist or scoring models in their R&D project funding process. 


The checklist approach involves the completion of profile charts for each project 
that is considered for funding. A list of criteria is first set to develop a profile. Then the 
candidate projects which meet these criteria or requirements are rated on a subjective scale 
such as high/medium/low or favorable/neutral/unfavorable. These ratings can be done by a 
committee or a single individual. The opinions of single individuals or committees can be 
summarized in a checklist by averaging their opinions. Figure 2.1 presents an example of a 
checklist. 

The advantage of a checklist lies in its simplicity and ease of use. A checklist can 
also provide a pictorial display of a project’s merits and limitations. Non-economical 
factors that are awkward or nearly impossible to include in more formal model 
constructions, such as social impacts and environmental concerns, can be added to a 


checklist without any difficulty. 


Criteria 
Profitability Marketability Success likelihood 
Poet A 


Project | 1” | sso all acl 
ee EEE 


Figure 2.1 Illustration of a Checklist For Three Projects [Ref. 4]. 





While the simplicity and the ease of the checklist model is very appealing, it can 
also be dangerous, since complex problems may be overlooked. Complicated relationships 
are not easily incorporated into checklists. Although many important factors may be 
included in a checklist, the relevance or importance of each individual factor is not 
captured, 

Scoring models attempt to remedy this problem by assigning weights to individual 
criteria and summarizing results as a single project score. Each candidate project is scored 
on each criterion, using an appropriate scoring scale (e.g., 10=excellent, 1=unacceptable). 
Each criterion is weighted relative to its importance. These scores and weights are 


combined according to the following model: 
= * 
fi, a » Ww ) Si 
j 


Here, 7; is the overall project score for project 1,w, is the relative criterion weight for the 
j" criterion, and Ss, is the criterion score for i" project on the j" criterion. Figure 2.2 


presents an example of a scoring model. 


Criterion Project Criterion 


Cniteria (1) a (wj) —_— (Si) Score ae Si) 


Probability of Success 
Profitability : 


Cost 2 
Patentability Rating ] 





Figure 2.2 Illustration of a Scoring Model [Ref. 4]. 


After finding the overall score for each project, the projects are prioritized from 
highest to lowest score. The highest scores are the most preferable to the decision maker. 
The weights and scores can be derived by using several methods. These include rank- 
ordering of attributes, pairwise comparisons, and other various rankings. Moore and 
Baker [Ref. 5] demonstrated that increasing the number of scoring intervals improves the 
accuracy of the scoring models. However, using psychometric testing, nine has been found 
as the maximum number of intervals that should be used [Ref. 5]. 

Moore and Baker [Ref. 6] compared scoring models with more sophisticated 
models such as economic, risk analysis, and constrained optimization models 
(mathematical programming models). Scoring models require less data input than the three 
other models. Because of data requirements, the scoring models have been viewed as most 
applicable during the initial stage of research, while the other model types are more 
appropriate for advanced research and engineering development analyses. They concluded 
that scoring models were not as constrained or weak as the three other models. 

Scoring models retain the advantage of checklists and profile charts in terms of 
their ability to consider a wide range of economic as well as non-economic criteria. 
Another advantage of the scoring model is the opportunity to use simple, low-cost 
methods for data acquisition. The scoring model also allows the decision maker to 
predetermine the impact of every factor while making his decision. Some examples of 
scoring model applications have been presented by Moore and Baker [Ref. 6], Dean and 
Nishry [Ref. 7], and Motley and Newton [Ref. 8]. 

The major disadvantage of the scoring model is that the project score is 
dimensionless which limits its use to rank-order comparisons. Another problem of scoring 
models is that the model development is non-formal. The weights and subjective scores 
can be determined in a number of ways, so there can be some difficulties in developing 


them precisely across different raters. 


B. FINANCIAL ANALYSIS MODELS 


Liberatore and Titus [Ref. 3] reported that almost every firm in their study was 
familiar with financial project selection and funding techniques and 62% of the firms they 
studied used these techniques. The two major methods of financial modeling were Net 
Present Value (NPV)/Internal Rate of Return (IRR) and benefit/cost analysis models. 

The IRR of a project can be defined as the effective rate of interest which will be 
earned on the money invested in the project. The NPV of a project is calculated by 
subtracting its present discounted cost from its present discounted benefit. A positive NPV 
shows that initial and future costs of project are less than the benefits of project, and it is 
desirable. The projects with the highest IRR or NPV are preferred. 

Benefit cost ratio models attempt to measure an estimated benefit or return from a 
project consistent with cost. Costs and benefit values are quantified in discounted dollar 
values. Costs are total resource costs of supporting the research project and benefits are 
the net earnings to be realized from the project. If the benefit/cost ratio of a project is 
greater than or equal to one, we can accept that project; otherwise, we can ignore It. 

Benefit cost ratios can be easily expanded to include a wide range of 
considerations such that a number of risk factors which reduce the expected project 
benefits. These risk factors can be shown as the probability of success of the project at 
various stages of development. Jackson [Ref. 9] presented an example of a benefit/cost 
ratio model, Olsen’s model, as: 

Td * nse ap an 
Total Project Cost 


where V is the economic value of the project, s is the annual sales volume derived from 
the project if the project succeeds, p is the profit per unit, and r, d, and m are the 
probabilities of research, development, and marketing success, respectively. The product’s 
expected life span is represented by n. Souder [Ref. 10] presented another index model, 


Ansoff’s model: 


r*d*m*(T+B)*E 


Fi fee — 
el Total Project Cost 


where r, d, and m are the same as Olsen’s model. T and B are subjective ratings of the 
technical and business merits of the project, respectively, and E is the present value of the 
earnings expected if the project succeeds. 

The above benefit/cost ratio models do not include other subjective (non- 
economic) benefits or costs, such as social, environmental or political effects. These 
benefits and costs can be entered into the model by expressing them in dollar values like 
the other factors. 

Benefit/cost ratio models are desirable since they overcome the dimensionless 
problem of checklist models and scoring models. There are, however, definite 
shortcomings of this model. One problem with this type of model is that the measures of 
expected values (inputs) are very difficult to obtain. The probabilities and cost and benefit 
estimations are hard to measure and require specific experience on the subject. It is also 
hard to express many non-economic factors in dollar terms. The benefits of military 
project selection are especially difficult to express in dollar terms. 

Another problem with benefit/cost ratio models is that the benefit cost ratios are 
not a useful tool for evaluating the consequences of alternative funding levels. Each 
element of benefits in the ratio must be reevaluated if the funding level is reduced or 


increased. Finally benefit/cost ratio models do not recognize resource constraints. 
C. RISK ASSESMENT MODELS 


One of the more useful tools used in R&D project selection and funding 1s the risk 
assessment technique. Liberatore and Titus [Ref. 3] found that the most of the 
respondents in their survey were familiar with risk assessment techniques, however only 
35% of the respondents used these techniques in their project selection or funding process. 
The decision tree models and Monte Carlo simulation models (risk analysis models) are 


the two most important and applicable risk assessment techniques. 


Decision tree methods attempt to project the chain of activities that will occur 
between the beginning and the completion of a project. Each of the steps along the chain 
contains estimated values of probabilities of success, costs, and returns. By this means it 1s 
possible to predict the expected profitability of alternative projects. 

The decision tree structure consists of decision points or nodes, outcome nodes, 
and branches emanating to and from each node. Each branch or path has a certain 
outcome and risk associated with it. The decision tree model is built by constructing a 
pay-off matrix that contains all outcomes and probabilities. The optimum path 1s found by 
Starting at the right-hand side of the tree and following an expected value algorithm 
folding back to the starting point. At each node the expected value is calculated for all the 
branches leaving that node, and the path with highest expected value is the one to be 
selected. 

One of the main advantages of the decision tree models is that they are analytically 
simple and can be graphically presented. The graphical representation of decision trees 
helps to clarify the available strategies and the potential risks, regrets, and trade-off. It 
makes the decision tree an excellent communication tool when communicating with high- 
level managers. 

The applicability of decision tree models in R&D project funding situations was 
demonstrated by Jackson [Ref. 9] and Raiffa [Ref. 11], and an example of decision tree 
models for project selection was presented by Flinn and Turban [Ref. 12]. The major 
disadvantage of decision trees is that the outcomes at each node are represented as a few 
discrete events rather than a continuous distribution of possible outcomes. 

Models that depend on inputs that are influenced by chance or estimated with 
uncertainty are called stochastic. The stochastic procedure for evaluating density functions 
(the probability distribution of a project’s outcome or the rate of return) has become 
known as the Monte Carlo simulation technique. Monte Carlo simulation was introduced 
by John von Neumann and Stanislaw Ulam when they both worked on the Manhattan 


Project at the Los Alomos National Laboratory. 
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Monte Carlo simulation is based on the decision tree model. The difference 
between these two models is that in Monte Carlo simulation each of the nodes is replaced 
with a probability distribution and this provides a stochastic decision tree. The probability 
distribution at each decision point will likely affect the project’s outcome. The simulation 
model selects one random value from each of these decision point distributions and 
computes the rate of return for this particular combination of random values. Numerous 
sets of random values can be selected to show the alternative rates of return for each set to 
the decision maker. 

The Monte Carlo simulation technique generally provides a more accurate 
description of the R&D decision process and offers a better basis for making project 
selection than other methods. The improvements in outcomes also bring a significant 
increase in information requirements. The benefits from the successful completion of the 
project must be assessed and the probability distributions for each unknown research 
project outcome must be estimated, and these processes are difficult and costly in most 
situations. Hespos and Strassman [Ref. 13] developed the most renowned application of 
the Monte Carlo simulation technique to a R&D project selection and funding problem. 

The major disadvantage of these two risk assessment models is that both models 
do not deal with resource constraints. Like financial analysis models, these two models fail 


to allocate scarce resources among various research projects. 
D. MATHEMATICAL PROGRAMMING MODELS 


Mathematical programming models have been used to solve many resource 
allocation, project selection, and capital budgeting problems for three decades. 
Surprisingly, the Liberatore and Titus study [Ref. 3] found that although most of the 
respondents have familiarity with mathematical programming models, there was no usage 
of these models for R&D project selection in the firms that they surveyed. The types of 
mathematical programming that have been used in R&D project selection and funding are 
linear programming, non-linear programming, integer programming, and_ goal 


programming. 


The linear programming techniques are concerned with the efficient use or 
allocation of limited resources among two or more activities or projects to meet desired 
objectives. Asher [Ref. 14] gives an example of a linear programming model in a R&D 
project funding problem. The general form of that example and other project funding 


models are shown below as: 


maximize en (2-1) 
subject to > AX <B (2 -2) 
O< X <1 (2 -3) 


where X is an n-component vector representing the funding levels of projects, C is an n- 
component vector representing the contribution (net profitability) of various projects, A is 
an m xX n matmx representing resource usage of the projects, and B is the particular 
resource levels. The range of X represents the funding levels of the projects and varies 
from zero to one. In many situations, projects are selected for full funding or they are not 
selected at all. In these situations, X can only take on values of zero or one. 

While many functional relationships in a mathematical model may be linear in 
nature, some relationships realistically are nonlinear. Nonlinear programming models are 
similar to linear models; the difference between the two models is that the objective 
function, constraint equations or both are nonlinear. Tyler, Moore and Clayton [Ref. 13] 
give an excellent example of nonlinear programming approaches to the project funding 
problem. 

Linear, integer, and nonlinear programming models are restricted to the 
establishment of only a single objective function. Most real world problems involve 
multiple goals. Goal programming, a modification and extension of linear programming, 
was introduced to solve this problem by Charnes and Cooper [Ref. 16] in 1961. 

The goal programming approach allows us to meet a system of complex objectives 
rather than a single objective. Unlike linear programming, the goal programming objective 
function usually does not contain choice variables. Instead it contains the positive and 


negative deviation variables from the designated goal constraints. The objective function 
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tries to minimize these deviations, based on the relative importance (weights) assigned to 
them. 


The formulation can be expressed as follows: 


minimize Wd +Wod™ (2 - 4) 
subjectto > gX-d* +d> =G (2 -5) 
a= B (2 -6) 

a ome = 0 (2 -7) 


The variables d@ and d represent the negative and positive deviations from the goal 


constraints. The goal constraint coefficients and goals are shown as g and G, respectively 


(equation 2-5). Wand W~ are nonnegative constants representing the relative weights to 
be assigned to the respective negative and positive deviation variables. Equations 2-6 and 
2-7 are the same as linear programming constraints regarding resource constraints, 
resource availability, and nonnegativity. 

The goal programming models are most appropriate for modeling multi-attribute 
decision problems. Most of the methodologies used to solve linear programming 
problems, like duality, sensitivity analysis, nonlinear programming, etc. will work on goal 
programming problems with minor revisions to the algorithms. Salvia and Ludwig [Ref. 
17] describe a goal programming model at the Lord Corporation to solve the project 
funding problem. Schniederjans [Ref. 18] listed goal programming applications by 


categorizing them according to their subjects. 
E. PREFERRED MODEL 


In this chapter, the four general types of project funding models were discussed: 
subjective models, financial analysis models, risk assessment models, and mathematical 
programming problems. Though each model type is important and has appropriate 
applications, the goal programming model is most applicable model for determining the 
force structure of the TNSG. There are several reasons for using this model. The main 
reason for using the goal programming model is that it permits specification of multiple 


objectives or targets to be achieved. Goal programming also has a great deal of flexibility 
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that is lacking in the other models mentioned in this chapter. Finally, the goal 
programming algorithm allows the use of resource constraints. Ruefli [Ref. 19] gives an 
example of the goal programming approach in the U.S. Department of Defense to solve a 


force structuring problem similar to the. TNSG problem. 
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Ii. MODEL DEVELOPMENT 


The goal programming technique discussed in the last chapter requires objective 
function coefficients. The intent of this chapter is to determine a methodology for 
converting subjective judgments into objective function coefficients (weights) for the goal 
programming model. These coefficients or weights represent the importance of each ship 
type for each goal or constraint in the goal programming model formulation. The problem 
the decision maker faces is how to select a method for determining meaningful weights 


among the numerous available approaches. 


A. TRADITIONAL METHODS OF QUANTIFYING SUBJECTIVE 
JUDGMENTS 


Data are produced using some system of measurement. Stevens [Ref. 20] defines 
measurement as “ the assignment of numbers to observations according to a set of rules.” 
The numerical values assigned when measuring objects implies a scale of measurement. 
Stevens [Ref. 20] gives the four levels of measurement represented by the four types of 
scales: nominal, ordinal, interval, and ratio. A nominal scale is simply categorizing a set of 
data into mutually exclusive subclasses. In an ordinal scale, the numbers are assigned to 
the various instances of the property, so that the order of numbers corresponds to the 
order of magnitude of the instances. Interval scales have all the properties of an ordinal 
scale, and in addition the distances between any two numbers on the scale are equally 
spaced. A ratio scale has all the characteristics of an interval scale, but it also has a true 
zero point as its origin and has the property of proportionality. 

These levels of measurement are used in the process of manipulating the numbers 
into a meaningful value of the object being measured. This value would then be used as an 
objective function coefficient for the goal programming model. The four main traditional 


methods of quantifying subjective judgments are presented below. 
1. Numerical Rating Method 


The numerical rating method is a very simple and direct method for quantifying 


subjective judgments. This method was first proposed by Stevens [Ref. 21] as a method of 
1 


obtaining comparative rankings in psychophysical experiments. Judges are asked to 
associate rated items with fixed reference points. This can be done by assigning numbers, 
or by plotting points on a continuous number line. Lodge [Ref. 22] gives an example of 
numerical rating that is used in obtaining opinions about how serious certain crimes are. 
The crime of a stolen car parked on the street was selected as a reference point and given 
a seriousness score of 100. Subjects were asked to rate other crimes using this reference 
point. If the judges think another crime is one tenth as less serious than the car theft, the 
seriousness score of that crime should be 10; if the crime was considered two times more 
serious, the score of that crime should be 200. After obtaining all the subjective 
responses, the numerical estimates are computed by taking their geometric mean. The 


continuous number line representation of this example 1s shown in Figure 3.1. 


less serious crime car theft more serious crime 





Figure 3.1 Numerical Rating Continuous Line. 
The primary advantage of the numeric rating method is computational simplicity. 
On the other hand, there are some disadvantages of this method; the main problem 1s that 
there is no natural origin for judgments, and judges frequently disagree with the reference 
point positioning. There are also no bounds on the interval scale; the lower bound (no 


crime situation) of the crime example is set to zero, but there is no upper bound. 
2. Categorical Judgment Method 


A commonly used means of obtaining numerical ratings from subjective ratings 1s 
the categorical judgment method, wherein judges assign instances to previously ordered 
categories. For example, program managers could be asked to rate a project milestone risk 
according to a scale of very low, low, average, high, or very high. The number of 
categories used can range from two to nine according to desires and resources of the 
surveyor and the skills of the raters [Ref. 23]. 

The procedure begins by rating the item in question and then constructing a matrix 


of cumulative frequency data with n row instances and m column categories. Each entry of 
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this matrix represents the number of judges who rated instance i in category j. A 
cumulative relative frequency matrix is developed from this matrix. In this matrix, entries 
are the proportions of the judges rating instance 1 in or below category j. The elements of 
this matrix are considered as areas under a standard normal curve and are converted to the 
corresponding z values. These values are then recorded in a z;; matrix consisting of n rows 
and m-1 columns, since the rightmost column may be omitted for computational purposes. 
The row average, r;, and column average, c;, are calculated, and a grand average G is 


found by averaging the column averages. A column sum of squares is computed as shown 


below: 
ssc=d(c,-G) (3-1) 
j 
Then for each row, the following 1s computed: 
SSR = (z,-1) | (3-2) 
; 
The scale value of the instances, s;, are found by solving the following equation in 
each row: 


s =G-r*J/SSC/SSR. (3-3) 


The main disadvantage of the categorical judgment method is_ that it obtains 
values with interval scale properties that can be linearly transformed to any other scale. 
Although it is more sophisticated than the numerical rating method, it 1s still 
computationally easy. The major disadvantage of this method is that the precision of its 


results is limited by the number of categories selected for use in the survey [Ref. 23]. 
3: Least Squares Method 


The least squares method was first proposed by Guilford [Ref. 24] as a means of 
obtaining scaled interval values from ordinal or comparative judgments. The inputs of this 
procedure are obtained by asking judges to do some ordinal ranking of various instances of 
a selected property. 

The method is started by taking the responses of the judges comparing several 
items with respect to a particular characteristic. For example, a group of editors might be 
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asked to rate three word processing programs of different software companies in terms of 
the overall quality of the program. Suppose that an editor judges program B is better than 
program C, which is better than program A. These responses can be shown in a frequency 


matrix (Table 3.1). 





Table 3.1 Least Squares Method Scoring Matnix. 


Since B is the preferred program by the judge, the entries are made in the 
corresponding rows of the column B that were rated inferior to B, which in this case are 
rows A and C. Since C was the second preferred program, an entry is made in the C 
column and A row. 

The responses of all judges are recorded in this manner, and collected to a fj 
frequency matrix as in Table 3.2. The sum of the cross-diagonal elements of this matrix will 
be equal to the total number of the judges. For example, there were 100 judges in Table 


3.3, and, as an example 54 judges A superior to C and 46 judges C superior to A 


f A B C 
A = 28 46 
B 72 : 65 
C 54 35 = 


Table 3.2 Least Squares Method Observed Frequency Matrix. 


The next step is to convert the frequency matrix fj, to a probability matrix Pi. The P 


matrix can be calculated by using the following equation: 


ie 
P, = > (3-4) 


fe 
For the above example, the probability matrix was obtained and is shown in Table 3.3. 
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0.28 


0.50 
0,35 
Table 3.3 Least Squares Method Probability Matrix. 
The important thing in Table 3.3 is that the diagonal entries of the probability 





matrix are set equal to 0.5. The probability matrix is then converted to the standard normal 
matrix X by subtracting the mean value of 0.5 from each value of p, and dividing this 
differences by the standard deviation of pj. The Xj values are the standard normal 
variables corresponding to the Pj; values of the probability matrix. In our example, the Xj; 
matrix is shown in Table 3.4. The least squares estimate of scale values s; was obtained by 
taking the mean of each column in the matrix X. The least squares solution requires that 
all elements of matrix X be present, however, some entries of this matrix can be vacant. 
When probability matrix entries are equal to 1.00 or 0.00, we can not obtain 
corresponding Xj entries. Torgerson [Ref. 25] stated that other least squares procedures 
have been developed for both this situation and the situation which judges do not rank all 


instances (incomplete probability matrices). 





Table 3.4 Least Squares Method Standard Normal Matrix. 


As in the categorical judgment method, the least squares estimate of scale values 
are linearly transformable to other scales. The use of an ordinal rating scale requires less 


time and effort on the part of the judges than other methods, so the survey forms of this 
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method are simple to use. The judges can simply list the instances in the order of 
importance regarding the compared factor. The main disadvantage of this method 1s that it 
requires a large number of judges to produce a reasonably accurate probability matrix. The 
judges sometimes do not rank all instances for various reasons, and this makes the least 


squares procedure difficult to use in scale development. 
4. Constant Sum Method 


The constant sum method, developed by Comrey in 1950 [Ref. 26], quantifies 
subjective ratings using pairwise comparisons. In this method, each instance 1s compared 
with each other by splitting 100 points. There will be n(n-1)/2 pairs. that must be 
considered, and 100 points will be divided between each in accordance with absolute ratio 
of the greater to the lesser. For example, if a judge gives 80 points to instance A and 20 
points to instance B, this indicates that A is four times more important than B. In the same 
manner a split of 60-40 would indicate a ratio of three to two, and 50-50 that two 
instances have the same magnitude. 

An example is used to illustrate this procedure. Suppose two judges are asked to 
evaluate three books on the basis of their content. Table 3.5 represents their respective 


comparison matrices where pj; 1s the number of points is given to book 1 when compared 


ce 


eo 





A | A 30 40 
B 80 50 60 B 70 50 ~~ 80 
C 70 40 #450 C 60 20 #50 


Table 3.5 Constant Sum Method Comparison Matnix. 


with book j. Both judges preferred book A to B, book A to C, and book C to B, but the 
intensities of these endorsements are different. The next step is to construct a matrix V by 


averaging the p; values across judges as shown in Table 3.6. 
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Table 3.6 Constant Sum Method Average 
Comparison Matrix. 
Another matrix W is formed from matrix Vi. The w, values are computed using 
the following equation: 
Wij — Vij / Vii (3-5) 


The W matrix for this example is shown in Table 3.7. 


1.00 0.33 0.54 
3.00 1.00 253 


1.86 0.43 1.00 
Table 3.7 Constant Sum Method W Matrix. 





The scale values can be computed by taking the n™ root of each column product, 


where n 1s the number of instances compared in our problem. In other words, 


lin 
S,= TI | ; (3-6) 
The calculation and results of our book example are demonstrated in Table 3.8. 
S,= [(1.00)*(3.00)*(1.86)]'” = 1.732 


S. = [(0.33)*(1.00)*(0.43)]'? = 0.521 


S; = [(0.54)*(2.33)*(1.00)]"* = 1.079 





Table 3.8 Constant Sum Method Scale Value Calculation. 


The constant sum method provides quantitative values that allows linear 
transformations and arithmetic operations. The scale, therefore, presses ratio properties 


rather than interval scale properties. Consistency is one of the problems of this method 
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that occurs when the number of instances is large. The ability of judges must be greater in 
order to provide consistent ratings. 

The traditional methods of quantifying subjective evaluation are not adequate for 
determining the force structure of the Turkish Naval Surface Group. The disadvantages of 
these methods indicate that another method that does not suffer from these disadvantages 


is necessary. 
B. ANALYTIC HIERARCHY PROCESS (AHP) 


The Analytic Hierarchy Process, a method for quantifying subjective variables, was 
developed and introduced by Professor Thomas L. Saaty of the Wharton School, 
University of Pennsylvania, and documented in his book of the same name [Ref. 27]. In 
the past few years AHP has been used in various decision-making models. Vargas [Ref. 
28] lists over 25 specific applications of AHP. 

The AHP procedure that has been popularized by Saaty contains four basic steps. 
These four steps in converting subjective judgments of a decision problem into numerical 


values are shown in Figure 3.2. 





AHP PROCEDURE 





Step 1: Break decision problem into | 


hierarchical levels 


Step 2: Collect input data by pairwise 
comparisons of decision elements 


Step 3: Using eigenvalue method to estimate 


the relative weights of decision elements 


Step 4: Aggregate the relative weights at 
each level 


Figure 3.2. AHP Procedure. 
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The first step is to break down the problem into a hierarchy of decision elements. 
The decision maker must develop a logical representation of the factors and levels 
involved in the problem scenario. The top level of the hierarchy, called the focus, consists 
of only one element: the overall objective of the problem. The subsequent levels of 
hierarchy include more specific objectives, attributes, or factors necessary to achieve the 
overall objective. Details of these attributes increase at the lower levels of the hierarchy. 
The final level of hierarchy contains the specific decision alternatives. 

The following car buying problem will be used to illustrate the AHP procedure in 
detail. For the decision to purchase a new automobile, the objective of buying a car is 
placed at the top of the hierarchy. Attributes of the car that influence the buying decision, 
such as cost, dependability, and comfort, are placed in the next level of the hierarchy, and 


the various car alternatives are listed at level three of Figure 3.3. 


Level | 


Level 2 


Level 3 





Figure 3.3 Car Buying Hierarchical Scheme. 


In the second step, judges are asked to make pairwise comparisons of the factors 
of one level that contribute to achieving the objective of the next higher level using the 


pairwise comparison scale shown in Table 3.9. 
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Intensity of 


Importance Definition Explanation 
] Equal importance. Two elements contribute equally to the 


property 
Moderate importance Experience and judgment slightly favor 
of one over another. one element over another 


Essential or strong Experience and judgment strongly favor 
importance one element over another 


Very strong An element 1s strongly favored and its 
importance dominance is demonstrated in practice 


Extreme importance The evidence favoring one element over 
another 1s of the highest possible order 
of information 


2,4,6,8 Intermediate values Compromise 1s needed between two 
judgments 


Reciprocals | When activity i compared to j 1s assigned one of the above numbers, 
then activity } compared to 1 1s assigned its reciprocal. 





Rationals Ratios arising from forcing consistency of judgments. 


Table 3.9 AHP Pairwise Companison Scale [Ref. 29]. 


As in the constant sum method, each judge must make n(n-1)/2 pairwise 
comparisons (where n is the number of elements on a level of the hierarchy). Because of 
this need, Saaty has recommended keeping the number of elements in any level at no more 
than nine and the number of levels between three and five [Ref. 30]. For the car buying 
example, a comparison matrix in Table 3.10 was constructed to compare the three cars 
with respect to comfort. The table shows that, when compared to the Honda, the Toyota 


is one half as comfortable, and the Ford is one fourth as comfortable. 


Toyota Honda 





Table 3.10 AHP Car Example Comparison Matnix. 
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In step three, the relative weight of each element in a level is computed using the 
eigenvalue solution technique. These weights are found by normalizing the pairwise 
comparison matrix, summing over the rows and obtaining an average row sum. The row 
sums are the values of the priority vector (eigenvector). The priority vector of comfort for 
the car buying example is (0.143',0.286,0.571)'. These values state that the Honda is 
considered more comfortable than the other two cars, since the Honda comfort factor is 
twice the size of the Toyota and four times as large as the Ford value. 

The “comfort” matrix used in the car buying example is consistent. In other 
words, the responses of the judges are not contradictory or conflicting (the Honda is 
preferred to the Toyota, the Toyota is preferred to the Ford, and the Honda is preferred to 
the Ford). An inconsistent matrix, for example, would show that while the Honda is 
preferred to the Toyota and the Toyota preferred to the Ford, the Ford would be 
preferred to the Honda. Complex decisions with more levels, attributes, and choices often 
contain inconsistencies. AHP measures the overall consistency of judgments by means of a 
consistency ratio (CR) that shows the quality of the data that has been input in the 
comparison matrices. The CR is computed by first finding the consistency index (CI). The 
consistency index is determined using the following equation: 
(a 
(a-1) 


where /_,, is the largest eigenvalue of the observed matrix of pairwise comparison and n 


CI = Ca 


is the number of the elements. 

The consistency ratio can then be found by dividing the CI by a random index (RI). 
The RI are average consistency indices of randomly generated matrices whose reciprocal 
entries were taken at random from the values 1/9, 1/8, ........ Pi ae ee , 8, 9. Saaty 
{Ref. 29] found that the value of the consistency ratio should be 10 percent or less. If it 1s 
more than 10 percent, the judgments are considered inconsistent. The consistency ratio for 


this example is zero. 


' For example 0.143 is obtained by summing first row elements and dividing this value by the sum of all 


matrix elements. 
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Step four of the AHP aggregates the relative weights of the various levels from the 
previous step in order to produce a vector of composite weights. This vector is the 
weighted rankings of decision alternatives (or selection choices) with respect to the factor 
being studied. Step four starts at the top of the hierarchy and determines the weights at 
that level. These weights are then multiplied by the eigenvector at the next lower level 
obtaining new vectors. Repeating this procedure yields relative weights of the elements at 
the lowest level of the hierarchy. 

Illustrating step four in our car buying example, the composite priority of the cars 
with respect to all the criteria is obtained by multiplying the priorities of the cars under 
each criterion by the priority of the criterion (in our example, the criteria in increasing 
order are: comfort, dependability, and cost?) and adding them across criteria. The 
calculation of the eigenvectors of cost and dependability are not shown. The procedure 
involved in step four 1s demonstrated in Table 3.11. 

Step four of the car buying example shows that the Honda should be purchased 
based on the buyer’s vehicle preferences regarding cost, dependability, and comfort. The 


Ford and the Toyota are essentially equivalent in preference. 


Level 2 eigenvalues: Cost =0.5 
Dependability =0.3 
Comfort = 0.2 


Level 3 eigenvalues 


Cost 0.400 
Dependability 0.100 
Comfort 0.286 


Ford weight : (0.5)(0.4)+(0.3)(0.2)+(0.143)(0.2) = 0.289 
Toyota weight : (0.5)(0.4)+(0.3)(0.1)+(0.286)(0.3) = 0.316 
Honda weight : (0.5)(0.2)+(0.3)(0.7)+(0.2)(0.571) = 0.424 





Table 3.11 Car Buying Example Step 4 Demonstration. 


* The priority of criterions can be obtained by constructing a pairwise comparison matrix over criterions 
and calculating the weights in the same way shown in comfort example. 
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AHP has several advantages over the traditional methods of quantifying subjective 
judgments discussed earlier in this chapter. AHP requires the decision maker to logically 
structure a complex problem. AHP gives scaled values that are on a ratio scale, and AHP 
is the only method that provides a mechanism for checking on the consistency of the input 
data which is a material requirement if meaningful results are desired. The only 
disadvantage of this method is that it is computationally more complex. But this 
disadvantage can be easily eliminated using special software packages prepared to solve 
AHP applications. 

There are also more sophisticated approaches in decision theory such as multi- 
attribute utility theory (MAUT). In recent years, several articles about the disadvantages 
of AHP and superiority of MAUT over AHP have been published. Dyer [Ref. 31] points 
out that AHP suffers rank reversal (an alternative that is chosen as the best over a set X is 
not chosen when some alternative, perhaps an unimportant one, is excluded from X). He 
concluded that this problem can be solved by changing ratio scales with interval scales as 
in MAUT. Perez [Ref. 32] gave the comparison of the two methods and stated: 


One would expect MAUT, since it requires only the construction of an interval 
scale, to be suitable for a wider range of applications than AHP. However, one would 
also expect that AHP, since it builds a ratio scale, would be more suitable to some 
situations in which the subjacent structure had a strong distributive component, 
particularly those in which the coefficients of the distribution were not strongly 
affected by changes in the set of available alternatives. 


The other shortcoming stated by Dyer [Ref. 31] is the scaling method of AHP. The 
replies to these criticisms by Saaty [Ref. 33], Harker and Vargas [Ref. 34], and the 
corresponding counterreplies show that no consensus has been reached. There are also 
some advantages of AHP over MAUT. Davies [Ref. 35] stated, “AHP is the only 
decision-making methodology dealing formally with inconsistency of judgments and 1s, 
therefore, superior to multi-attribute utility.” 

Lee and Ahn [Ref. 36] used the AHP method for selecting ground warfare 
weapons for the Korean Army similar to the TNSG problem. AHP also has been 


successfully applied to many resource allocation problems similar to the problem of 


Jy) 


determining the force structure of TNSG. It is the method selected here for quantifying 


subjective judgments. 
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IV. GOAL PROGRAMMING MODELS 


There are three basic steps required in formulating a goal program. These steps 
are: define the decision variables and constants; formulate the constraints (system 
constraints and goal constraints); and develop the achievement (objective) function. This 
chapter will describe two main mixed-integer, goal programming (GP) models developed 
using these three steps. 

Model 1 is formulated to determine the number and mix of ship types to purchase 
or build given a constant budget and weighted AHP preferences for each ship type. In this 
model, the main goal is to maximize the effectiveness of the force mix when a constant 
budget is given. Model 2’s formulation uses the weapon requirements needed to meet the 
threat’s specific capabilities (defined in the mission areas of air, surface, submarine, and 
amphibious operations), ship weapon capacities for these requirements, and the weighted 
AHP preferences for the ship types. In the second model, the main goal is to minimize the 
cost of providing a fixed force effectiveness. After presenting the details of the models, we 


describe the motivation in this case for the three steps in model formulation stated above. 
A. DETAILED MODEL FORMULATIONS 


The goal programming models 1 and 2 are shown in Figures 4.1 and 4.2 below by 
using the above methodology. In subsequent chapters, model 1 will be referred to as the 


budget constrained, AHP-preferred (or budget/AHP) model. 
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Achievement Function 

Minimize: 3 (AHP Priority Weight for Ship Type 1 * 
Negative Deviation for Ship Type 1) 
= Total Deviation 

Goal Constraint 

Subject to: 

X; - Positive Deviation from Number of Ships of Type 1 Required 
+ Negative Deviation from Number of Ships of Type 1 Required 
= Unattainable Goal 

System Constraint 


>», (X; *Cost;) < Total Budget 


X, Positive and Negative Deviations are Integer Variables > 0 





Figure 4.1 Model 1 Goal Programming Formulation. 
In NPS format, Model 1 ts as follows: 
Indices : 
i: Type of combat ships (CG, DDG, FFG, FSG, PBFG); 
w: Priority weights; 
Data: 
COST; : Cost of each ships of type1 in billions of dollars; 
WEIGHTSiw : AHP weights for each ship type; 
NUM; . Desired number of ships of type 1 (unattainable); 
TOTBUDGET : Total budget available; 
Decision Variables : 
X; . Number of ships of type 1 to purchase; 
DEVNEG; : Negative deviation from desired number of ships of type 1; 
DEVPOS; : Positive deviation from desired number of ships of type 1; 
DEVIATION : Deviation from objective function; 


Model Formulation : 


Minimize DEVIATION = >| DEVNEG, * > WEIGHTS,, (Objective function) 


Subject to X, - DEVNEG, + DEVPOS = NUM, Vi (Goal constraints) 


> COST, * X, < TOTBUDGET (System constraints) 
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X,, DEVNEG,, DEVPOS, = 0. and integer 


Figure 4.2 presents Model 2. In subsequent chapters, it will be referred to as the 
weapon-constrained model. The AHP priority weight for ship type is summed across the 


AHP weight for all categories. 
Achievement Function 


Minimize: a (AHP Priority Weight for Ship Type 1 * 
Negative Deviation for Ship Type 1) 
= Total Deviation 
Goal Constraint 
Subject to: 
Xj - Positive Deviation from Number of Ships of Type i Required 
+ Negative Deviation from Number of Ships of Type 1 Required 
= Unattainable Goal 


System Constraint 

» (X; *capacity for Long Range SSM) = Lower Limit of Long Range SSM Required 

> (X; *capacity for Long Range SSM) < Upper Limit of Long Range SSM Required 

> (X; *capacity for Short Range SSM) = Lower Limit of Short Range SSM Required 
Ss (X; *capacity for Short Range SSM) < Upper Limit of Short Range SSM Required 
DS (X; *capacity for Long Range SAM) = Lower Limit of Long Range SAM Required 
yy (X; *capacity for Long Range SAM) < Upper Limit of Long Range SAM Required 
y (X; *capacity for Short Range SAM) = Lower Limit of Short Range SAM Required 
yy (X; *capacity for Short Range SAM) < Upper Limit of Short Range SAM Required 
Dy (X; *capacity for Long Range ASM) = Lower Limit of Long Range ASM Required 
(X; *capacity for Long Range ASM) < Upper Limit of Long Range ASM Required 
(X; *capacity for Short Range ASM) = Lower Limit of Short Range ASM Required 
yy (Xj; *capacity for Short Range ASM) < Upper Limit of Short Range ASM Required 
(X; *capacity for Long Range Guns) = Lower Limit of Long Range Guns Required 
a (X; *capacity for Long Range Guns) < Upper Limit of Long Range Guns Required 

> (Xj; *capacity for Short Range Guns) = Lower Limit of Short Range Guns Required 
x 


(X; *capacity for Short Range Guns) < Upper Limit of Short Range Guns Required 
Positive and Negative Deviations are Integer Variables 2 0 


Figure 4.2 Model 2 Goal Programming Formulation. 
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In NPS format, Model 2 is as follows: 

Indices : 

i: Type of combat ships (CG, DDG, FFG, FSG, PBFG); 

p: Parameters for each ship type (weapon types which each ship type contains); 
w: Pniority weights; 

Data : 

VALUES; : parameter values for each ship type; 

WEIGHTS, : AHP weights for each ship type; 

NUM; . Desired number of ships of type 1 (unattainable); 

LOWREQ, : Lower limit (number) of launchers required from weapon type p; 
UPREQ, : Upper limit (number) of launchers required from weapon type p; 
TOTBUDGET : Total budget available; 

Decision Variables : 

X; . Number of ships of type 1 to purchase; 
DEVNEG; . Negative deviation from desired number of ships of type i; 
DEVPOS; : Positive deviation from desired number of ships of type 1; 
DEVIATION : Deviation from objective function; 


Model Formulation : 


Minimize DEVIATION = >. DEVNEG, * >| WEIGHTS,, (Objective function) 


Subject to X,- DEVNEG, + DEVPOS = NUM, Vi (Goal constraints) 
>_ VALUES, *X, 2 LOWREQ, Vp (System constraints) 
>_VALUES,,*X,L<UPREQ,  p (System constraints) 


X,, DEVNEG,, DEVPOS, 20 and integer 


B. DECISION VARIABLES AND CONSTANTS 


Determination of the decision variables and might hand side is the first step in 


construction of a goal programming model. The nght hand side constants may be either 
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resource constraints or specified goal levels. Decision variables are controlled by the 
decision maker and are sometimes referred to as “control variables”. Decision variables 
also must be nonnegative. In both models, the main goal is to find the optimum number 
and the optimum mix of ship types. In these GP models the decision variables, denoted by 
X(i), represent the numbers of ships of type i to purchase (i= CG, DDG, FFG, FSG, 
PBFA). These vanables are nonnegative integer variables. 

In each goal constraint, there is a negative and positive deviation variable that 
represents the amounts by which the goals are overachieved (positive deviation) or 
underachieved (negative deviation). These are also nonnegative integer variables. 

In Model 1, constants are the numbers of the ship types desired and the total 
budget. The total budget 1s the amount which will be used to structure a new Naval 
surface force in the next five years. The numbers of ships desired were set as unattainable 
goals for both models. Model | ensures that the total budget is spent in order to reach the 
goal. Model 2 ensures that all weapon requirements will be met. These weapon 
requirements are the numbers of weapon launchers for specific missions to meet the threat 


countries’ capabilities. 
C. CONSTRAINTS 


When analyzing the relationships among decision variables and their relationships 
to the goals, a set of constraints should be formulated. There are two types of constraints 
in a goal programming formulation, system constraints and goal constraints. System 
constraints are absolute constraints, and the decision maker has no control over these 
constraints. System constraints must be satisfied before an optimal or satisfactory solution 
can be considered. 

Goal constraints are not absolute constraints. These constraints contain positive 
and negative deviation variables, because goal programming attempts to minimize the 
deviation from the goal rather than attempt to satisfy the goal completely. The solution 
could result in overachieving or underachieving the particular goal. 

Each goal constraint in Model | requires an unattainably high number of each ship 


type. This will result in a negative deviation from the desired (unattainable) amounts in 


a8 


order to force a positive deviation which can be minimized in the objective function. The 
objective function will try to minimize this positive deviation. The deviation is the penalty 
for underachieving the desired number of ships which is weighted by the ship’s AHP 
preference weight. Model 2 has exactly the same goal constraints as Model 1. 

The system constraint in Model 1] 1s concerned with the total money available (our 
budget) for the ship types selected. This constraint states that the total cost of the selected 
ship mix can not exceed our total budget. The total budget 1s used to structure a new 
Naval surface force over the next five years. The system constraints in Model 2 are 
different from the Model 1. These constraints require specific numbers of weapon 
launchers for each mission (AAW, ASUW, ASW, and Amphibious Warfare, also divided 
into short and long range). Values for these requirement levels will be determined by the 
Turkish Naval Surface Group Staff by examining the force structure and weapon 
capabilities of the threat countries. But in this formulation, dummy values are used as 
requirement values and budget. Hence, this model is a generic model and can be used by 
any country by changing these requirements and budgets. 

A goal programming model requires that all decision variables be greater than or 
equal to zero. The negative and positive deviation variables and all other variables used in 


the goal constraints must be nonnegative. 
D. OBJECTIVE FUNCTION 


The final step in the model development is to state the objective function. All 
constraints and goals must be completely identified in the model in order to analyze each 
goal in terms of whether over or underachievement of the goal is satisfactory or not. If 
overachievement is satisfactory (as in Model 1 and 2), the positive deviation variables can 
be eliminated from our objective function. Likewise, if underachievement is acceptable, the 
negative deviation variables can be eliminated from objective function. If exact 
achievement is required, then both negative and positive deviation variables will be in the 
objective function. In both models the objective function minimizes the sum of negative 


deviations from the number of ships of varying types weighted by the AHP priority vector. 
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E. VARIATIONS ON BASIC MODELS 


A variation of Model 1 can be formulated to represent a least cost, war situation 
(minimum deviation from the budget). The scenario on which weapon requirements are 
based is a conventional war situation with a specific country. This model formulation can 
be used to determine the least cost mix of ship types to meet the total weapon 
requirements in a war situation. Values for weapon requirement levels will also be 
determined as in Model 2. This model will not use any AHP weights, and it is simply a 
linear programming formulation. 

The objective functions of the models used are linear (first-order) functions. But 
linear utility functions can give irrelevant results when there are strong interactions 
between factors (this interdependency problem is solved in our models by using structural 
constraints). Quadratic (second-order) objective function models can also be used to solve 
this interaction problem. The AHP models are modified using a quadratic objective 


function and are evaluated in our determination of the force mix of TNSG. 
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V. SOFTWARE SUPPORT AND MODEL RESULTS 


The primary objective of this study is to develop a model and to use this model to 
determine the optimum force structure of TNSG. This chapter is intended to present the 
computer programs used to solve the model and the model results. Specifically, this entails 
a discussion of the data collection and software development process, the numerical 


output acquired from the various models’ iterations and the outcomes produced. 
A. DATA COLLECTION AND SOFTWARE DEVELOPMENT 


The process of collecting data and designing software to support the GP model 
formulation is an essential step in determining force structure of the TNSG. As discussed 
earlier, the GP model requires that the objective function weights be determined by the 
AHP. The AHP survey was intended to determine a relative measure of capability of the 
ships in several warfare situations and factors. The survey was based on recommendations 
by Saaty [Ref. 29] for eliciting preferences of paired comparisons. The forms were 
distributed to the 20 Turkish Navy Officers attending the Naval Postgraduate School. 
Fifteen of the surveys were completed. It 1s assumed that those officers who did complete 
the survey are qualified to make these judgments. A copy of the actual surveys can be 
found in Appendix A. 

Computer software was used to perform the necessary AHP calculations and 
determine the objective function coefficients of the GP model. The Microsoft Excel V.7.0 
spreadsheet and S-Plus V.3.3 Release 1 statistical package by Mathsoft [Ref. 37] were 
chosen to perform AHP calculations. Excel especially, is easy to use in many statistical 
calculations. The eigenvalues were calculated by Excel. The eigenvectors and consistency 
ratios were calculated by S-Plus. 

The Generalized Algebraic Modeling System (GAMS 2.25) [Ref. 38] was selected 
to solve goal programming models. GAMS is designed to ease the construction and 
solution of large and complex mathematical programming models. The programs are 


formulated as mixed integer problems (MIP) and require the use of the XA subsolver of 
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GAMS that is the implementation of the Sunset Technology XA Callable Library. XA 
contains high performance solvers for LP and MIP. The XA solver first solves the problem 
as a linear program. Then it uses the Pivot and Complement heuristic to find an initial 
integer feasible solution. Finally XA, uses a branch and bound algorithm to determine 


whether there is a better solution and to verify optimality [Ref. 39]. 
B. MODEL RESULTS 


The output from the AHP procedure and the GP models will be surveyed in this 
section. This evaluation will include a comparison of different GP models and selection 


of an optimal model. 
1. AHP Results 


The AHP surveys were collected and the pairwise comparison matrices were 
obtained by taking the geometric mean as recommended by Saaty [Ref. 29] . These 
matrices were put on Excel and S-Plus to generate the final AHP results. The equal 
preference weights are given for all ship characteristics in computation of the AHP 
results. The AHP comparison matrices for each ship characteristic (cost, speed, warfare 
capabilities, surveillance capability, and fuel consumption) and AHP results (AHP 2 


values, consistency indices, consistency ratios and the eigenvectors) are presented in 


Appendix B. Table 5.1 shows the composite matrix of the AHP values. 


cf . Re AAW | ASW | ASUW | AMPW | CIWS | Fuel Surv. 
ip ost pee Cap Cap. Cap. Cap. Cap. Con. Cap. 





Table 5.1 Composite Matrix of AHP Results. 


Consistency ratios for all comparison matrices were well under the recommended 


value of 10 percent. The highest consistency ratio was only 3.5 percent. It is critical to the 
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GP model that the standard for consistency not be violated, since inconsistency will 
adversely affect the accuracy and credibility of the model. The CG was preferred in five of 


nine categories. 
2: Goal Programming Results 


As indicated earlier, dummy data values are used in this study to make these 
models generic and usable for other countries. Appendices C through G show the GAMS 
formulation and summarized output. 

Appendix C contains the formulation and results of Model 1. It is a budget- 
constrained, AHP-preferred mix goal to determine the mix and the number of ship types. 
The resulting mix consisted of 20 PBFA and 1 FSG; no other types of ships are included. 
All of the budget is used; there was no remaining budget. When the AHP preferences are 
compared with the ship costs, the PBFA is the most preferred ship type in the model. A 
variation of Model 1 was formulated to represent a war situation, least cost (minimum 
deviation from the budget) goal. The model formulation is used to determine the least cost 
mix of ship types to meet the total weapon requirements in a war situation. This model is 
shown in Appendix D, and the ships selected were 2 DDGs, 6 FFGs, 4 FSGs, and 1 
PBFA. The multi-mission capability and cost (when compared with the other multi- 
mission ship types) of the FFGs are the main reasons for the selection of this type. 

Model 2 (Appendix E) was formulated using a weapon constrained, AHP- 
preferred mix goal. The AHP preferences and weapon requirements drove the selected 
mix of ships. This mix was 8 FFGs, 2 FSGs, and 3 PBFAs. The addition of AHP 
preferences increased the number of FFGs and PBFAs. 

Table 5.2 summarizes the model results. As can be seen, the objective value of the 
budget/AHP model is the smaller than the AHP-preferred model (the same goals and 
objective function). The appearance of all PBFAs in the budget/AHP solution is explained 
by the high ratio of AHP weight values to the cost for this ship type. This results in the 
minimum deviation value in this model. But PBFAs are not multi-mission surface 


combatants, and the mix with only this type of ship seems inappropriate for this problem. 
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Weapon-Constrained Models 


Ship Type | Budget AHP | Minimum Cost AHP Preferred 
CG os 





Table 5.2 Gams GP Model Results Summary. 


The evaluation of the weapon-constrained models can be done by comparing the 


ship types on the basis of cost and benefit. Table 5.3 summarizes the results of this 


comparison. The sum of each ship type’s eigenvalues are normalized to one 


Marginal Benefit 


Total Mix 


Total Mix 
Cost ($107) 


Total Mix 
Benefit 


Benefit-Cost 
Ratio 





Table 5.3. Weapon-Constrained Models Comparison. 


and used as a measure of marginal benefit. It can be seen that the AHP-Preferred mix, 
while costing 0.8 billion dollars less, has approximately equal total benefit with the 
minimum cost model and has a larger benefit-cost ratio. The minimum cost model is the 


only formulation not using AHP weights, and it has the least benefit-cost ratio. Our 
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evaluation of the models shows that the AHP-Preferred model gives the most appropriate 
mix for the TNSG. These results provide an indication of the overall superiority of FFG- 
dominated ship mixes. The cost-benefit ratio calculations were made using results of the 
AHP-Preferred and weapon constrained model. The use of the AHP cost weights in the 
cost-benefit calculations may be viewed as away of accounting for cost when firm cost 
estimates are not available. Future research on this issue is meredit. 

The budget-constrained AHP model can be modified using a quadratic objective 
function with the same hypothetical AHP values. The GAMS formulation of this new 
model and its output are shown in Appendix F. Unlike Model 1 (budget-constarined, 
AHP-preferred), this quadratic programming (QP) model was not dominated by PBFAs 
(Model 1 mix contains 20 PBFA and 1 FSG, QP model mix contains 4 FFG, 4 FSG, and 5 
PBFA). QP results are more similar to weapon-constrained model and they can be used in 


future studies. 
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VI. SENSITIVITY ANALYSIS 


In previous chapters, three different models were developed based on several data 
values and objective function weights obtained from the surveys. It is obvious that 
responses to these surveys may not be perfect; the data used may be subject to error, and 
ship capabilities can change with time. The purpose of this chapter is to analyze the 
impacts of changes in the goal programming model. There will be sensitivity analyses of 
the following changes: total budget, weapon requirements, and objective function 


coefficients (weights). 
A. CHANGES IN BUDGET 


Model 1 (budget-constrained and AHP-preferred) is designed to find a surface 
force mix with a budget constraint and AHP weights. The possible budget for force 
structuring 1s a five year plan, and it can change over time due to several effects (political, 
economic, etc.). Possible budget levels were run in the model and resulted in different 
solutions. The solutions showed that the number of ships increase (decrease) with an 
increase (decrease) in the budget, but the ship types that are selected by the model did not 
change. 

The least-cost model was also run at several different budget levels. Every budget 
level gave the same solution to the model even with a changing objective function. The 
objective function value, the weighted deviation, increases (decreases) with decreases 
(increases) in the budget. It showed that the least cost model is not sensitive to changes in 


the budget. 
B. CHANGES IN WEAPON REQUIREMENTS 


Weapon requirement constraints force specific numbers of weapon launchers for 
each mission (AAW, ASUW, ASW, and Amphibious Warfare, also divided into short and 
long range). Values for these requirement levels will be determined by the Turkish Naval 
Surface Group Staff after examining the force structure and weapon capabilities of the 


threat countries. These requirements can change with a force structure change of the 
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threat countries and/or addition of new threats (for example, the threat country can decide 
to reduce or increase its force size). 

The minimum weapon requirements must first be satisfied for all warfare areas. We 
will try to determine the effects of changes in the weapon requirements. Every constraint 
value is changed using a loop in the GAMS model. The following results were found. 

The analysis of the requirements for long and short range surface to surface 
missiles (NLSMUP: the upper bound for the number of long range surface to surface 
missiles, NLSMLO: the lower bound for the number of long range surface to surface 
missiles, NSSMUP: the upper bound for the number of short range surface to surface 
missiles, and NSSMLO: the lower bound for the number of short range surface to surface 
missiles) gave the results that are shown in Table 6.1. Table 6.1 shows that changes in the 
lower bounds of these constraints do not have any effects on the solution of the problem. 
These constraints are likely to be redundant. However, changes in the upper bounds have 


significant effects on the determination of the force mix. Decreases inthe upper bounds 


Right Hand Side Obj. Func. Value Force Mix 


157 22 8 FFG, 2 FSG, 3 PBFA 
160.48 1 DDG, 6 FFG, 3 FSG, 1 PBFA 
NO FEASIBLE INTEGER SOLUTION 


ld ie22 8 FFG, 2 FSG, 3 PBFA 
NO FEASIBLE INTEGER SOLUTION 


157.22 8 FFG, 2 FSG, 3 PBFA 
158.23 1 DDG, 7 FFG, 2 FSG, 2 PBFA 
159.85 2 DDG, 6 FFG, 2 FSG, 1 PBFA 
161.80 2 DDG, 5 FFG, 3 FSG 

NO FEASIBLE INTEGER SOLUTION 


157.22 8 FFG, 2 FSG, 3 PBFA 
NO FEASIBLE INTEGER SOLUTION 





Table 6.1 Sensitivity Analysis of Long and Short Range Surface to Surface Missile 
Requirements. 
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result in lower numbers of ships in the force mix which results in bigger objective function 
values. 

The analysis of the requirements for long and short range surface to air missiles 
(constraints: NLAMUP, NLAMLO, NSAMUP, and NSAMLO) gave the results that are 
shown in Table 6.2. Table 6.2 shows that a change in the upper bound of long range 
surface to air missile (SAM) requirements does not have any effect on the solution of the 
problem. Again, this constraint is likely to be redundant. The changes in the lower bound 
of long range SAM and upper bound of short range SAM requirements have significant 
effects on the determination of the force mix. The increase in the lower bound of short 
range SAM does not have any effect in the force mix, but at some point the problem 


becomes infeasible. 


Right Hand Side Obj. Func. Value Force Mix 


NLAMUP re. 22 8 FFG, 2 FSG, 3 PBFA 
NO FEASIBLE INTEGER SOLUTION 


NLAMLO ise 8 FFG, 2 FSG, 3 PBFA 
158.58 1 DDG, 7 FFG, 1 FSG, 3 PBFA 
159.85 2 DDG, 6 FFG, 2 FSG, 1 PBFA 
NO FEASIBLE INTEGER SOLUTION 


a 22 8 FFG, 2 FSG, 3 PBFA 
boron 1 DDG, 7 FFG, 2 FSG, 2 PBFA 
158.58 1 DDG, 7 FFG, 1 FSG, 3 PBFA 
NO FEASIBLE INTEGER SOLUTION 


157.22 8 FFG, 2 FSG, 3 PBFA 
NO FEASIBLE INTEGER SOLUTION 





Table 6.2 Sensitivity Analysis of Long and Short Range Surface to Air Missile 
Requirements. 


Table 6.3 shows the analysis of the changes in the long and short range 
antisubmarine weapon requirements (constraints: NLUASMUP, NLASMLO, NSASMUP, 
and NSASMLO). A change in the upper bound of long range antisubmarine weapon 
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requirements does not have any effect on the solution of the problem; this constraint 1s 
also likely to be redundant. The changes in the other antisubmarine weapon requirements 
have significant effects on determination of the force mix. The lower bound of short range 


antisubmarine weapons results in very sensitive mix changes. 


Right Hand Side | Obj. Func. Value Force Mix 


li a2? 8 FFG, 2 FSG, 3 PBFA 
NO FEASIBLE INTEGER SOLUTION 


157.22 8 FFG, 2 FSG, 3 PBFA 
158.54 1 DDG, 7 FFG, 2 FSG, 2 PBFA 
159.58 1 DDG, 7 FFG, 1 FSG, 3 PBFA 
158.54 1 DDG, 7 FFG, 2 FSG, 2 PBFA 
NO FEASIBLE INTEGER SOLUTION 


hie 2 8 FFG, 2 FSG, 3 PBFA 
158.58 1 DDG, 7 FFG, 1 FSG, 3 PBFA 
NO FEASIBLE INTEGER SOLUTION 


NSASMLO See 6 bRGeZ PSGno Bien 
se) 8) 1 DDG, 6 FFG, 3 FSG, 2 PBFA 
NO FEASIBLE INTEGER SOLUTION 





Table 6.3 Sensitivity Analysis of Long and Short Range Anti Submarine Weapons 
Requirements. 

Table 6.4 shows the analysis of the changes in the requirements of long and short 
range guns (constraints: NUGUNSUP, NLGUNSLO, NSGUNSUP, and NSGUNSLO). 
Table 6.4 shows that changes in the upper bounds of long and short range guns have no 
effects on the solution of the problem. These constraints are also likely to be redundant. 
The changes in the lower bounds of long and short range guns have significant effects on 


the force mix. These requirements are likely to be dominant constraints that affect the 


number of other weapons. 
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Right Hand Side | Obj. Func. Value Force Mix 


NLGUNSUP 22 8 FFG, 2 FSG, 3 PBFA 


NLGUNSLO 157.26 8 FFG, 1 FSG, 4 PBFA 
157,22 8 FFG, 2 FSG, 3 PBFA 
157.26 8 FFG, 2 FSG, 4 PBFA 
SOD 8 FFG, 2 FSG, 3 PBFA 
158.54 1DDG,7FFG,2 FSG, 2 PBFA 
159.85  2DDG, 6 FFG, 2 FSG, 1 PBFA 
NO FEASIBLE INTEGER SOLUTION 


NSGUNSUP 722 8 FFG, 2 FSG, 3 PBFA 


NSGUNSLO 156.4 10 FFG, 3 PBFA 
156.58 9 FFG, 2 FSG, 2 PBFA 
157,22 8 FFG, 2 FSG, 3 PBFA 
159.13 1 DDG, 6 FFG, 3 FSG, 2 PBFA 
NO FEASIBLE INTEGER SOLUTION 





Table 6.4 Sensitivity Analysis of Long and Short Range Guns Requirements. 

As a summary, the sensitivity analysis of weapon requirements show that the 
changes in the lower bounds of long and short range guns give the most sensitive results. 
Decreases in the upper bound of the requirements decreases the number of the ships in the 
force mix. The increases in the lower bound of the long range guns forces the model to 
include the ship types CG and DDG (the ship types have more long range guns), and this 
increases the number of long range SSMs (because CG and DDG have more long range 


SSM). This makes the changes in the upper bound of the long range SSM more sensitive. 


GC: CHANGES IN OBJECTIVE FUNCTION COEFFICIENTS (WEIGHTS) 


Survey results were converted, using the AHP, into objective function coefficients. 
The survey is based on subjective judgments that can possibly vary for several reasons. In 
this section an analysis of the impact of the changes in the subjective evaluations will be 
conducted. 

First, a sensitivity analysis can be done to determine the impact of a mistake in 


completing a survey or transcribing data from a survey. This analysis is illustrated in Table 
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6.5, a table that displays an actual matrix taken from the AHP output in Appendix B and a 
flawed matrix that contains an input error. The only change that was made in the flawed 
matrix is that the upper nght-hand number was changed from 0.17 to 1.0, representative 
of a common input error. The consistency ratio (CR) computation shows that this small 
mistake increases the CR to 0.168. Since a CR this high (bigger than 0.1) is unacceptable, 
the matrix data input would have to be examined, and the error would be corrected. This 


example shows the sensitivity of the GP model in responding to minor lapses. 


Accurate Matrix Flawed Matrix 
[OO3 0500250 200017 00, 0/5050'255 0:20 MiEeo 
200" T0080 50 0-25 70-20 . §2:300 TGOR0 50 90 25020 
400 2.00 1.00 0.50 0.33 4.00 2.00 1.00 0.50 0.33 
5.00 4.00 2.00 1.00 0.51 5.00 4.00 2.00 1.00 0.51 
5.95 5.00 2.99 1.95 1.00 1.00 5.00 2.99 1.95 1.00 
Consistency Ratio: 0.024 Consistency Ratio: 0.168 


Table 6.5 Comparison of Accurate Matrix with Flawed Matrix. 


The GP model might also be subjected to a change of opinion. A respondent to a 
survey could decide that a ship type was judged inappropriately. This can be shown using 
an example of such a situation. Consider the comparison matrix of ship types according to 


cost as shown in Table 6.6. 


CG DDG 
100 0.50 
Z.00. iiree 
400 2.00 
5.00 4.00 
>> 5.00 


Consistency Ratio: 0.024 
Eigenvector: .049 .091 





Table 6.6 Omnginal Comparison Matrix for Cost. 
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Suppose the surface group staff officers receive information implying that new 
DDGs will have both gas turbine and diesel engines that will increase their cost. The 


officers might then change the original comparison matrix to one shown in Table 6.7. 


CG DDG 
CG EOC) O75 
DDG Fso8 100 
FFG 4.00 2.50 
FSG 5.00 4.44 
Piiawes 5) 5.50 


Consistency Ratio: 0.012 
Eigenvector: .054 071 





Table 6.7 Comparison Matrix for Cost with Minor Changes. 


The only changes were made in the second row and second column of the new 
matrix, values corresponding to the DDG. The numbers in the orginal matrix were larger, 
signifying the cost of DDGs is not too great; the lower numbers in the new matrix denote 
that the cost of DDGs has increased significantly. This is also reflected by the coefficient 
eigenvector. The largest coefficient change occurred in the DDG value, which dropped 
from 0.091 to 0.071, but the other coefficient weights increased only slightly. 

The example above demonstrates the sensitivity of the models to major subjective 
changes in the responses to the AHP surveys. The new coefficient weights were used in 
the AHP-preferred models, and the results show that big changes in subjective judgments 
make big differences in the surface group force mix. These changes have even greater 


effects in the quadratic GP models. 
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VU. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


The determination of force structure of a naval surface group is essentially a 
project selection problem. There are additional factors, aside from cost and weapon 
requirements, which must be included when considering the mix and size of the naval 
surface group. These factors should include not only the ship attributes such as speed, 
cost, warfare capabilities, and fuel consumption, but with respect to these attributes, the 
relative advantage held by each ship in accomplishing the mission. A survey is a method to 
obtain subjective judgments rating ships with respect to the above factors. The Analytic 
Hierarchy Process appears to be a useful tool for converting these subjective judgments 
into numerical preference weights. 

Goal programming formulations using AHP-weighted deviation variables can 
produce a consistent, and in some sense, the best solution for determining force size and 
mix. This was demonstrated by the smallest goal deviation and the highest benefit-cost 
ratio achieved by the ship mix selected by the AHP-Preferred model formulation. 

Although CGs and DDGs are the most capable ships, the mix does not contain 
these ship types because of their high cost. FFGs are the most preferable ship types, 
because they have multi-mission capability and low cost when compared with other multi- 
mission ship types (CG and DDG). FSGs and PBFAs are preferable because of their low 


cost, but their warfare capabilities limit their numbers in the force mix. 
B. RECOMMENDATIONS FOR FUTURE STUDY 


Additional model formulations should be constructed using other surface warfare 
requirements, such as a helicopter landing capability and operational capability in open sea 
versus shoreline (i.e, seas with lots of islands like the Aegean Sea). This model 
formulation can be enlarged with the addition of the other warfare types, such as mine 


wartare. 
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The AHP survey should be expanded to include other factors which are thought to 
contribute to naval surface group mission effectiveness (in several warfare types). Some 
factors can be divided into subgroups (ASUW capability in open sea or shoreline). It 
might also prove useful to ascertain the relative importance of mission effectiveness 
parameters among themselves (e.g., the advantage of AAW warfare capability over fuel 
consumption). The mixed integer GP models can be converted to quadratic mixed GP 
models using a different applications of the AHP results. Additionally, the population of 
survey judges should be increased. 

Future studies can include more thorough analysis of AHP versus MAUT. Also, 
new methods can be obtained by combining the attractive features of both techniques 
(e.g., multiple judgments and calculation of inconsistency in AHP and the 
interdependencies in MAUT). 

The participation of the ship types in a task force group is another future issue in 
our study. A nominal task force must be a task-organized group of combatants and 
includes the interdependencies between ship types that are included in the force. The task 
group formation must provide all kind of warfare capability in the protection of the all 


constituent variables. 


a2 


APPENDIX A. SURVEY 
SHIP CHARACTERISTICS SURVEY 

This survey was designed to obtain pairwise comparison data of the ships in the 
naval surface group. The ships are compared on the basis of nine different characteristics. 
You have been briefed on the purpose of this survey and its theoretical foundations. 

The following pairwise comparison scale will be used to compare each type of ship 
against the others with respect to the characteristics: 
1 - equal/same 
2 
3 - moderate 
4 
5 - strong 
6 
7 - very strong 
8 
9 - extremely strong 


The values 2, 4, 6, and 8 are intermediate values. 


Please ensure that you circle either advantage or disadvantage when making each 
comparison. 


Example: The CG givesa__3 _advantage/disadvantage over the FFG. 


Questions may be directed to LTJG Erol Unal, 642-9860. 
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Comparison 1: Costs. Given the following costs (billions) for each ship type: CG - 


The CG gives a: advantage/disadvantage over DDG. 
advantage/disadvantage over FFG. 
advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


The DDG gives a: advantage/disadvantage over FFG. 
advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


The FFG gives a: advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


The FSG gives a: advantage/disadvantage over PBFA. 


Comparison 2: Speed. The speed of a ship is a factor that increases the operational 
capability of the ship in different warfare situations. In addition, speed is a factor in the 
survivability of the ship. 


The CG gives a: advantage/disadvantage over DDG. 
advantage/disadvantage over FFG. 
advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


The DDG gives a: advantage/disadvantage over FFG. 
advantage/disadvantage over FSG. 


advantage/disadvantage over PBFA. 


The FFG gives a: advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


The FSG gives a: advantage/disadvantage over PBFA. 

Comparison 3,4,5,6,7; Warfare Capabilities. This shows the capability of the ships in 
different warfare areas. The warfare capability of the ship changes according to the ship’s 
weapons in that warfare area. | 


With respect to AAW capability: (AAW weapons and countermeasures) 


The CG gives a: advantage/disadvantage over DDG. 
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advantage/disadvantage over FFG. 
advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


The DDG gives a: advantage/disadvantage over FFG. 
advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


The FFG gives a: advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


The FSG gives a: advantage/disadvantage over PBFA 
With respect to ASUW capability: (ASUW weapons and countermeasures) 


The CG gives a: advantage/disadvantage over DDG. 
advantage/disadvantage over FFG. 
advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


The DDG gives a: advantage/disadvantage over FFG. 
advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


The FFG gives a: advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


The FSG gives a: advantage/disadvantage over PBFA. 
With respect to ASW capability: (ASW weapons and countermeasures) 


The CG gives a: advantage/disadvantage over DDG. 
advantage/disadvantage over FFG. 
advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


The DDG gives a: advantage/disadvantage over FFG. 
advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


The FFG gives a: advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


The FSG gives a: advantage/disadvantage over PBFA. 


a) 


With respect to AMPW capability: (Navy gun fire support capability) 


The CG gives a: advantage/disadvantage over DDG. 
advantage/disadvantage over FFG. 
advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


The DDG gives a: advantage/disadvantage over FFG. 
advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


The FFG gives a: advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


The FSG gives a: advantage/disadvantage over PBFA. 


With respect to CIWS capability: 


The CG gives a: advantage/disadvantage over DDG. 
advantage/disadvantage over FFG. 
advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


The DDG gives a: advantage/disadvantage over FFG. 
advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


The FFG gives a: advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


The FSG gives a: advantage/disadvantage over PBFA. 

Comparison 8: Surveillance capability. (Radars, ESM, and ECM capability) 

With respect to Surveillance capability: 

The CG gives a: advantage/disadvantage over DDG. 
advantage/disadvantage over FFG. 
advantage/disadvantage over FSG. 


advantage/disadvantage over PBFA. 


The DDG gives a: advantage/disadvantage over FFG. 
advantage/disadvantage over FSG. 
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The FFG gives a: 


The FSG gives a: 


advantage/disadvantage over PBFA. 


advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


advantage/disadvantage over PBFA 


Comparison 9: Fuel consumption. 


With respect to fuel consumption: 


The CG gives a: 


The DDG gives a: 


The FFG gives a: 


The FSG gives a: 


advantage/disadvantage over DDG. 
advantage/disadvantage over FFG. 
advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


advantage/disadvantage over FFG. 
advantage/disadvantage over FSG. 


advantage/disadvantage over PBFA. 


advantage/disadvantage over FSG. 
advantage/disadvantage over PBFA. 


advantage/disadvantage over PBFA 


i, 





CG 
DDG 
FFG 
FSG 
PBFA 


CG 
1.000 
2.000 
4.000 
5.000 
Seb 


LAMBDA : 4.892 


CONSISTENCY INDEX : 0.027 
CONSISTENCY RATIO : 0.024 


DDG 
0.500 
1.000 
2.000 
4.000 
5.000 


WEIGHT EIGENVECTOR : 


CG 
DDG 
FFG 
FSG 
PBFA 


CG 
1.000 
0.638 
0.519 
0.187 
Oa) 


LAMBDA : 4.942 


CONSISTENCY INDEX : 0.0145 
CONSISTENCY RATIO : 0.013 


DDG 
1.567 
1.000 
0.667 
0.250 
0.167 


WEIGHT EIGENVECTOR : 


CG 
DDG 
FFG 
FSG 
PBFA 


CG 
1.000 
O25 
0.704 
0.400 
0.333 


DDG 
1.380 
1.000 
0.899 
0.667 
0.476 


FFG 
0.250 
0.500 
1.000 
2.000 
2.994 


CG 
0.049 


FFG 
25 
1.500 
1.000 
O26! 
0.200 


CG 
O37) 


FFG 
1.420 
Peer 
1.000 
0.625 
0.455 


FSG 
0.200 
0.250 
0.500 
1.000 
1.953 


DDG 
0.091 


FSG 
5.340 
4.000 
34750 
1.000 
0.400 


DDG 
0.272 


FSG 
2.500 
1.500 
1.600 
1.000 
0.813 
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COMPARISON MATRIX FOR COSTS 


PBFA 
0.168 
0.200 
0.334 
Ole 
1.000 


FFG 
0.181 


PBFA 
8.240 
6.000 
5.000 
2.500 
1.000 


FFG 
e727 


PBFA 
3.000 
2.100 
2.200 
1.230 
1.000 


FSG 
0.289 


COMPARISON MATRIX FOR AAW CAPABILITY 


FSG 
0.087 


COMPARISON MATRIX FOR ASUW CAPABILITY 


APPENDIX B. AHP RESULTS FOR SHIP CHARACTERISTICS 


PBFA 
0.390 


PBFA 
0.039 


LAMBDA : 5.004 
CONSISTENCY INDEX : 0.001 
CONSISTENCY RATIO : 0.00089 


CG DDG FFG FSG PBFA 
WEIGHT EIGENVECTOR : vee O21 0.220 onss 0.106 


COMPARISON MATRIX FOR ASW CAPABILITY 


CG DDG FFG FSG PBFA 
CG 1.000 Os» 0.806 3.760 5.600 
DDG Eo 1.000 1.142 4.850 8.205 
FFG 1.24] 0.876 1.000 4.560 7.145 
BSG 0.266 0.206 wie 1.000 3.250 
PBFA 0.179 O22 0.140 0.308 1.000 
LAMBDA : 5.063 
CONSISTENCY INDEX : 0.0158 
CONSISTENCY RATIO : 0.0140 

CG DDG FFG FSG PBFA 

WEIGHT EIGENVECTOR : 0.238 O33] 0.297 0.099 0.035 


CG 
DDG 
FFG 
FSG 
PBFA 


CG 
1.000 
0.495 
0.340 
One? 
0.141 


LAMBDA : 5.082 


CONSISTENCY INDEX : 0.0205 
CONSISTENCY RATIO : 0.0183 


DDG 
2.020 
1.000 
0.649 
0267 
0.168 


WEIGHT EIGENVECTOR : 


FFG 
2.940 
1.540 
1.000 
32 | 
0.209 


CG 
0.388 


FSG 
5210 
37750 
3.120 
1.000 
0.398 


DDG 
0.270 
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PBFA 
7.100 
5.950 
4.780 
ZO 10 
1.000 


FFG 
0.210 


COMPARISON MATRIX FOR AMPW CAPABILITY 


RSG PBFA 
0.091 0.041 


CG 
DDG 
FFG 
FSG 
PBFA 


CG 
1.000 
0.890 
0.881 
0.475 
0.170 


LAMBDA : 5.001 


CONSISTENCY INDEX : 0.00025 
CONSISTENCY RATIO : 0.00022 


DDG 
iel2s 
1.000 
0.995 
0.528 
0.180 


WEIGHT EIGENVECTOR : 


CG 
DDG 
EEG 
FSG 
PBFA 


CG 
1.000 
0.508 
0.447 
0.266 
0.183 


LAMBDA : 5.062 


CONSISTENCY INDEX : 0.0155 
CONSISTENCY RATIO : 0.0139 


DDG 
IED G7 
1.000 
0.899 
0.320 
0.215 


WEIGHT EIGENVECTOR : 


CG 
DDG 
FFG 
EsG 
PBFA 


CG 
1.000 
1.124 
1.143 
1.845 
ZAS5 


DDG 
0.890 
1.000 
1.109 
1.321 
2.294 


FFG 
135 
1.005 
1.000 
0.580 
0.184 


CG 
0.287 


FFG 
2236 
bi 
1.000 
0.358 
0.236 


CG 
0.356 


FFG 
0.875 
0.902 
1.000 
als 
1.764 


FSG 
2.105 
1.895 
1.723 
1.000 
0.310 


DDG 
0.263 


FSG 
59/56 
D123 
Dt 
1.000 
0.448 


DDG 
OZe7 


FSG 
0.542 
Oat 
0.867 
1.000 
eon 
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PBFA 
5.890 
So o7 
5.425 
3230 
1.000 


FFG 
O25 


PBFA 
5.456 
4.654 
4.234 
2.234 
1.000 


EEG 
G 232 


COMPARISON MATRIX FOR SPEED 


PBFA 
0.363 
0.436 
0.567 
0.789 
1.000 


COMPARISON MATRIX FOR CIWS CAPABILITY 


FSG 
0.148 


COMPARISON SURVEILLANCE CAPABILITY 


FSG 
0.103 


PBFA 
0.047 


PBFA 
0.051 


LAMBDA : 5.018 
CONSISTENCY INDEX : 0.0045 
CONSISTENCY RATIO : 0.0040 


CG DDG FFG FSG PBFA 
WEIGHT EIGENVECTOR : OBZ Cay 0.169 0.220 rons 


COMPARISON MATRIX FOR FUEL CONSUMPTION 


CG DDG FFG FSG PBFA 
CG 1.000 0.512 0.418 0.398 0.332 
DDG 1.355 1.000 0.884 0.689 0.496 
FFG 22 inst 1.000 0.857 0.568 
FSG 213 1.45] Pel67 1.000 0.612 
PBFA 3.012 2.016 1.761 1.634 1.000 
LAMBDA : 5.014 
CONSISTENCY INDEX : 0.0035 
CONSISTENCY RATIO : 0.0031 

CG DDG FFG FSG PBFA 

WEIGHT EIGENVECTOR : 0.089 0.169 0.200 0.226 0.316 
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APPENDIX C. GAMS GOAL PROGRAMMING FORMULATION AND 
OUTPUT FOR BUDGET-CONSTRAINED AND AHP PREFERRED MIX 


MODEL 


A. GAMS FORMULATION 


SETS 


I type of combat ships /CG cruisers 


DDG guided missile destroyers 

FFG guided missile frigates 

FSG guided missile corvettes 

PBFA guided missile fast patrol boats / 


P parameters for each ship type 
/C_ cost of ship type I in billions of dollars / 


W priority weights 


/CW 
SW 


AAWW 
ASWW 
ASUWW 
AMPW 


CIWS 
FW 
SCW 


TABLE 


cost weights 

speed weights 

aaw capability weights 

asw capability weights 

asuw capability weights 

amphibious warfare capability weights 
ciws capability weights 

fuel consumption weights 
surveillance capability weights /; 


VALUES (1,P) parameter values for each ship type 
C LSM SSM LAM SAM LASM LSSM LGUNS GUNS 


CG 
DDG 
FFG 
FSG 
PBFA 


TABLE 


ey 16 0 8 0 8 6 Z 2 
6 8 0 8 0 8 3 2 2 
“ 4 + 4 4 0 6 l 0 
20 4 4 0 4 0 6 2 
ne 4 +} 0 0 0 0 0 Z 


WEIGHTS (],W) AHP weights for each ship type 


CG 
DDG 
FFG 
FSG 


CW 


049 
O91 
181 
289 


SW AAW ASWW ASUWW AMPW CIWS FW SCW 


neo 6315 238 onl 388 287 089 .356 
els 2 22 331 2) wae, OSE OOme 25). 
HOS e277) 97] eo Ue 5S pee eae 
20 eO So! 1092 ioe 1091 148 .226 .103 
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PBEA 390 3727) G32 055 .106 04] 047 BiG 051 


PARAMETER 
NUM (I) — desired number of ship type I (unattainable) 
LCG 20 
DDG 20 
FFG = 20 
ESGwe20 
PBFA 20 /; 


SCALARS 
TOTBUDGET /2.8 /,; 


INTEGER VARIABLES 
X (I) number of ships I to purchase 
DEVNEG (I) negative deviation from desired number of ships of type I 
DEVPOS (I) positive deviation from desired number of ships of typeI ; 


POSITIVE VARIABLE 
DEVIATION deviation from objective function ; 

EQUATIONS 
BUDGET stay within budget 
NSHIPS(J) desired number of ship type I (unattainable) 
OBJECFUN objective function; 
BUDGET... SUM(, X(D * VALUES(I,’C’)) =L= TOTBUDGET ; 
NSHIPS()).. X(1) - DEVPOS(I) + DEVNEG(I) =E= NUM(D) ; 
OBJECFUN.. DEVNEG( ‘CG’ ) * SUM(W, WEIGHTS(‘CG’,W)) 


+ DEVNEG( ‘DDG?’ ) * SUM(W,WEIGHTS(‘DDG’,W)) 
+ DEVNEG( ‘FFG’ ) * SUM(W,WEIGHTS(‘FEG’,W)) 

+ DEVNEG( ‘FSG’ ) * SUM(W,WEIGHTS(‘FSG’,W)) 

+ DEVNEG( ‘PBFA’ ) * SUM(W,WEIGHTS(‘PBFA’,W)) 
=E= DEVIATION : 


MODEL GP /ALL/; 
SOLVE GP USING MIP MINIMIZING DEVIATION; 


B. GAMS MIXED INTEGER PROBLEM OUTPUT 


SOLVE SUMMARY 


MODEL GP OBJECTIVE DEVIATION 
TYPE > MiP DIRECTION MINIMIZE 
SOLVER XA FROM LINE 78 
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**** SOLVER STATUS 1 NORMAL COMPLETION 


**** NIODEL STATUS 1 OPTIMAL 
**** ORJECTIVE VALUE 151.5620 
RESOURCE USAGE, LIMIT 0.000 50000.000 
ITERATION COUNT, LIMIT 5 20000 
No better solution than: 151.56200 

Absolute Relative 
Actual distance 0.00000 0.00000 


Tolerances (OPTCA) 0.00000 (OPTCR) 0.00000 


LOWER LEVEL UPPER MARGINAL 


---- EQU BUDGET -INF 2.800 2.800 -6.990 

LOWER LEVEL UPPER MARGINAL 
---- EQU OBJECFUN é ; -1.000 
---- VAR X number of ships I to purchase 


LOWER LEVEL UPPER MARGINAL 


CG . 100.000 4.058 
DDG 100.000 2.168 
FFG 100.000 0.106 
FSG 1.000 100.000 
PBFA 20.000 100.000 


---- VAR DEVNEG negative deviation from desired number of ships of type I 


Orvek ) LEVEL UPPER MARGINAL 


CG . 20.000 100.000 

DDG : 20.000 100.000 

FFG ; 20.000 100.000 

FSG 19.000 100.000 
PBFA : 100.000 0.043 


---- VAR DEVPOS positive deviation from desired number of ships of type I 


LOWER LEVEL UPPER MARGINAL 
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CS : 100.000 Zee SS 


DDG . . 100.000 2076 
FFG , 100.000 oe 
PSG 100.000 1.398 
PBFA ; 100.000 0.909 


LOWER LEVEL UPPER MARGINAL 


---- VAR DEVIATION -INF 151.562 +INF 
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APPENDIX D. GAMS GOAL PROGRAMMING FORMULATION AND 
OUTPUT FOR THREAT-CONSTRAINED (LEAST COST) MODEL 


A GAMS FORMULATION 


SETS 


I type of combat ships / CG cruisers 


DDG guided missile destroyers 

FFG guided missile frigates 

FSG guided missile corvettes 

PBFA guided missile fast patrol boats / 


P parameters for each ship type 


ec 
LSM 
SSM 
LAM 
SAM 
LASM 
SASM 


cost of ship type I in billions of dollars / 

number of the long range surface to surface missiles 

number of the short range surface to surface missiles 

number of the long range surface to air missiles 

number of the short range surface to air missiles 

number of the long range anti submarine rockets or torpedoes 
number of the short range anti submarine rockets or torpedoes 


LGUNS number of the long range guns 
SGUNS number of the short range guns /; 


TABLE 


VALUES (I,P) parameter values for each ship type 
C LSM SSM LAM SAM LASM LSSM LGUNS GUNS 


CG 
DDG 
FFG 
FSG 
PBFA 


3 


SCALARS 


i’ 16 0 8 0 8 6 2 
6 8 0 8 0 8 3 Z 
3 4 Zz 4 4 0 6 ] 
2 + 4 0 4 0 6 ] 
al a 4 0 0 0 0 0 


TOTBUDGET (eos 1; 


POSITIVE VARIABLES 
DEVIATION deviation from objective function 
DEVNEGB negative deviation from budget 


Z 
2 
0 
2 
Z 


DEVPOSB positive deviation from budget; 
INTEGER VARIABLE 
X(I) number of ships I to purchase ; 
EQUATIONS 
NUMLSM number of Ism launchers needed to meet threat’s surface power 
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NUMSSM number of ssm launchers needed to meet threat’s surface power 


NUMLAM number of lam launchers needed to meet threat’s air power 
NUMSAM number of sam launchers needed to meet threat’s air power 
NUMLASM number of lasm launchers needed to meet threat’s submarines 
NUMSASM number of sasm launchers needed to meet threat’s submarines 


NUMLGUNS _ number of lguns needed to meet threat’s power 
NUMSGUNS _ number of sguns needed to meet threat’s power 


BUDGET stay within budget 

OBJFUN objective function; 

NUMLSM.. SUM(I, X(I) * VALUES(,’LSM’)) =G= 40 ; 
NUMSSM.. SUMC( I, X(I) * VALUES(,’SSM’)) =G= 40 ; 
NUMLAM.. SUMC(C I, X(1) * VALUES(I,’LAM”’)) =G= 40° ; 
NUMSAM.. SUMC( I, X(1) * VALUES(I,’SAM”’)) =G= 40 ; 


NUMLASM.. SUM(L, X(I) * VALUES(I,"LASM’)) =G= 30 ; 
NUMSASM.. — SUM(I, X(I) * VALUES(I,’SASM’)) =G= 30 ; 
NUMLGUNS.. SUM(I, X(I) * VALUES(I,,LGUNS’)) =G= 13 ; 
NUMSGUNS.. SUM(I, X(I) * VALUES(I,’SGUNS’)) =G= 13 ; 


BUDGET. SUM( I, X(I) * VALUES(I,’C’)) - DEVPOSB + DEVNEGB 
=E=TOTBUDGET: 
OBJFUN.. DEVPOSB =E= DEVIATION; 
MODEL GP /ALL/: 


SOLVE GP USING MIP MINIMIZING DEVIATION; 


B. GAMS MIXED INTEGER PROBLEM OUTPUT 


SOLVE SUMMARY 


MODEL GP OBJECTIVE DEVIATION 
TYPE MIP DIRECTION MINIMIZE 
SOLVER XA FROM LINE 74 


**** SOLVER STATUS 1 NORMAL COMPLETION 
**** MODEL STATUS 1 OPTIMAL 
=T** OBJECTIVE VALUE 1.1300 


RESOURCE USAGE, LIMIT 0.160 50000.000 
ITERATION COUNT, LIMIT 4 20000 


No better solution than: 1.13000 


Absolute Relative 
Actual distance 0.00000 0.00000 
Tolerances (OPTCA) 0.00000 (OPTCR) 0.00000 
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---- EQU NUMLSM 
---- EQU NUMSSM 
---- EQU NUMLAM 
---- EQU NUMSAM 
---- EQU NUMLASM 
---- EQU NUMSASM 
---- EQU NUMLGUNS 
---- EQU NUMSGUNS 
---- EQU BUDGET 
---- EQU OBJFUN 


---- VAR DEVIATION 
---- VAR DEVNEGB 


LOWER 


40.000 
40.000 
40.000 
40.000 
30.000 
30.000 
13.000 

13.000 

2.800 


LOWER LEVEL 


-INF 


BEVEL 


60.000 

44.000 
40.000 
40.000 
34.000 
48.000 
14.000 

14.000 

2.800 


50 


UPPER MARGINAL 


+INF 

+INF , 

+INF 07025 

+INF 0.050 

+INF 

+INF 

+INF 

+INF , 

2.800 -1.000 
-1.000 


UPPER MARGINAL 


+INF 
+INF 1.000 


---- VAR DEVPOSB 


ee WHO 


CG 
DDG 
Fata 
FSG 
PBFG 


REVEL 


2.000 
6.000 
4.000 
1.000 


number of ships I to purchase 


+INF 


UPPER MARGINAL 


100.000 
100.000 
100.000 
100.000 
100.000 


0.700 
0.400 


0.130 
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APPENDIX E. GAMS GOAL PROGRAMMING FORMULATION AND 
OUTPUT FOR THREAT CONSTRAINED AND AHP PREFERRED MODEL 


A. GAMS FORMULATION 


Ses 
I type of combat ships / CG 
DDG 
FFG 
FSG 
PBFA 


P parameters for each ship type 


cruisers 
guided missile destroyers 
guided missile frigates 
guided missile corvettes 
guided missile fast patrol boats / 


/C cost of ship type I in billions of dollars / 
LSM number of the long range surface to surface missile launchers 
SSM number of the short range surface to surface missiles launchers 
LAM _ number of the long range surface to air missiles launchers 
SAM number of the short range surface to air missiles launchers 
LASM number of the long range anti submarine rockets or torpedoes 
SASM number of the short range anti submarine rockets or torpedoes 
LGUNS number of the long range guns ( for aw and asuw ) 
SGUNS number of the short range guns ( for asuw and aaw ) / 


W priority weights 


/CW cost weights 
SW speed weights 


AAWW _ aaw capability weights 

ASWW _. asw capability weights 

ASUWW -asuw capability weights 

AMPW amphibious warfare capability weights 
CIWS ciws capability weights 

FW fuel consumption weights 

SCW surveillance capability weights ie 


TABLE 


VALUES (1,P) parameter values for each ship type 
C LSM SSM LAM SAM LASM LSSM LGUNS GUNS 


CG y 16 0 
DDG 6 8 0 
REG =D of “: 
FSG Zz 4 4 
PBFA .13 4 + 


8 0 8 6 y 2 
8 0 8 3 2 Z 
4 4 0 6 ] 0 
0 + 0) 6 ] 2 
0 0 0 0 0 2 


71 


TABLE 


WEIGHTS (I,P) parameter values for each ship type 


CW SW AAW ASWW ASUWW AMPW CIWS FW SCW 
CG 049 .132 .375 238 319 388 287 ©6089 356 
DDG “OS! -152> 2272 foul 22) 270 263 169 237 
FFG 188 21699 4227 297 220 a2 110 255 200 32382 
FSG .289 .220) =.087 .099 135 091 148 226 .103 
PBFA .390 327 .039 035 106 04] 047 316 .051 
PARAMETER 
NUM (I) | desired number of ship type I (unattainable) 
[CG 20 
DDG 20 
FFG 20 
FSG 20 
PBFA 20 /; 
INTEGER VARIABLES 
MD) number of ships I to purchase 
DEVNEG (1) negative deviation from desired number of ships of type I 
DEVPOS (I) positive deviation from desired number of ships of type I 
FREE VARIABLE 
DEVIATION _ deviation from objective function ; 
EQUATIONS 
NLSMUP number of Ism launchers needed to meet threat’s surface power 
(upper limit) 
NLSMLO (lower limit) 
NSSMUP number of ssm launchers needed to meet threat’s surface 
(upper limit) 
NSSMLO (lower limit) 
NLAMUP number of lam launchers needed to meet threat’s air power 
(upper limit) 
NLAMLO (lower limit) 
NSAMUP number of sam launchers needed to meet threat’s air power 
(upper limit) 
NSAMLO (lower limit) 
NLASMUP number of lasm launchers needed to meet threat’s submarines 
(upper limit) 
NLASMLO _ (lower limit) 
NSASMUP number of sasm launchers needed to meet threat’s submarines 
(upper limit) 
NSASMLO (lower limit) 


a2 


NLGUNSUP 
NLGUNSLO 
NSGUNSUP 
NSGUNSLO 
NSHIPS(1) 
OBJFUN 


NLSMUBP.. 
NLSMLO.. 
NSSMUBP.. 
NSSMLO.. 
NLAMUP.. 
NLAMLO.. 
NSAMUP.. 
NSAMLO.. 
NLASMUP.. 
NLASMLO.. 
NSASMUB.. 
NSASMLO.. 
NLGUNSUP.. 
NLGUNSLO.. 
NSGUNSUP.. 
NSGUNSLO.. 
NSHIPS(I).. 
OBJECFUN.. 


number of lguns needed to meet threat’s power (upper limit) 
(lower limit) 


number of sguns needed to meet threat’s power (upper limit) 
(lower limit) 

desired number of ship type I (unattainable) 

objective function; 


SUM( I, X(I) * VALUES(I,’LSM’)) =L= 55 
SUM( I, X(I) * VALUES(I,’LSM’)) =G= 30 
SUM( I, X(1) * VALUES(I,’SSM’)) =L= 55 
SUM( I, X(I) * VALUES(I,’SSM’)) =G= 30 
SUM(I, X(I) * VALUES(I,"LAM’)) =L= 55. ; 
SUM( I, X(1) * VALUES(I,’"LAM’)) =G= 30 ; 
SUM( I, X(I) * VALUES(I,’SAM’)) =L= 55 ; 
SUM( I, X(I) * VALUES(I,’SAM’)) =G= 30; 
SUM( I, X(I) * VALUES(I,’"LASM’)) =L= 41 

SUM( I, X(I) * VALUES(I,"LASM’)) =G= 23 

SUM( I, X(I) * VALUES(L,’SASM’)) =L= 41 

SUM( I, X(I) * VALUES(L,’SASM’)) =G= 23 
SUM(I, X(I) * VALUES(I,"LGUNS’)) =L= 15 
SUM( I, X(I) * VALUES(,’ LGUNS’)) =G= 11 
SUM(I, X(I) * VALUES(I,’SGUNS’)) =L= 15 
SUM(I, X(I) * VALUES(I,’SGUNS’)) =G= 11 
X(I) - DEVPOS(I) + DEVNEG(I) =E= NUM()) ; 
DEVNEG( ‘CG’ ) * SUM(W, WEIGHTS(‘CG’,W)) 


2 
2 


2 


> 
2 
2 


> 


+ DEVNEG( ‘DDG’ ) * SUM(W, WEIGHTS(‘DDG’,W)) 
+ DEVNEG( ‘FFG’ ) * SUM(W, WEIGHTS(‘FFG’,W)) 

+ DEVNEG( ‘FSG’ ) * SUM(W, WEIGHTS(‘FSG’,W)) 

+ DEVNEG( ‘PBFA’ ) * SUM(W, WEIGHTS(‘PBFA’,W)) 
=E= DEVIATION ; 


MODEL GP /ALL/,; 


2 


SOLVE GP USING MIP MINIMIZING DEVIATION; 


B. GAMS MIXED INTEGER PROBLEM OUTPUT 


SOLV 


MODEL GP 
TYPE = MiP 
SOLVER XA 


E SUMMARY 


OBJECTIVE DEVIATION 
DIRECTION MINIMIZE 
FROM LINE 121 


**** SOLVER STATUS — 1 NORMAL COMPLETION 


**** MODEL STATUS 


1 OPTIMAL 
nS 


**** OBIECTIVE VALUE 


RESOURCE USAGE, LIMIT 
ITERATION COUNT, LIMIT 


No better solution than : 


Actual distance 
Tolerances (OPTCA) 


---- EQU NLSMUP 
---- EQU NLSMLO 
---- EQU NSSMUP 
---- EQU NSSMLO 
---- EQU NLAMUP 
---- EQU NLAMLO 
---- EQU NSAMUP 
---- EQU NSAMLO 
---- EQU NLASMUP 
---- EQU NLASMLO 
---- EQU NSASMUP 
---- EQU NSASMLO 
---- EQU NLGUNSUP 
---- EQU NLGUNSLO 
---- EQU NSGUNSUP 
---- EQU NSGUNSLO 


157.22000 


Absolute 


0.00000 


lo72200 


0.330 
BS 


0.00000 (OPTCR) 


LOWER LEVEL UPPER MARGINAL 


-INF 
30.000 
-INF 
30.000 
-INF 
30.000 
-INF 
30.000 
-INF 
23.000 
-INF 
23.000 
-INF 
10.000 
-INF 
10.000 


52.000 

52.000 
52.000 
52.000 
32.000 
32.000 
40.000 
40.000 
24.000 
24.000 
36.000 
36.000 
10.000 
10.000 
10.000 
10.000 


50000.000 
20000 


Relative 
0.00000 
0.00000 


55.000 
+INF 
55.000 
+INF 
55.000 
+INF 
55.000 
+INF 
41.000 
+INF 
41.000 
+INF 
15.000 
+INF 
15.000 
+INF 


LOWER LEVEL UPPER MARGINAL 


---- EQU OBJFUN -1.000 
---- VAR X number of ships I to purchase 
LOWER LEVEL UPPER MARGINAL 
CG 100.000 -2.233 
DDG 100.000 -2.026 
FFG 8.000 100.000 EPS 
FSG 2.000 100.000 -1.398 
PBFA 3.000 100.000 -1,.352 
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---- VAR DEVNEG 


LOWER 


CG 
DDG 
FFG 
FSG 
PBFA 


---- VAR DEVPOS 


LOWER 


CG 
DDG 
FFG 
Sy 
PBFA 


---- VAR DEVIATION 


negative deviation from desired number of ships of type I 


EW EL 


20.000 
20.000 
12.000 
18.000 
17.000 


UPPER MARGINAL 


100.000 
100.000 
100.000 
100.000 
100.000 


leo 


positive deviation from desired number of ships of type I 


LEVEL UPPER MARGINAL 


100.000 
100.000 
100.000 
100.000 
100.000 


26 
2.026 


1.398 
1 352 


LOWER LEVEL UPPER MARGINAL 


-INF 


e220 


Ws 


+INF 
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APPENDIX F. GAMS QUADRATIC GOAL PROGRAMMING FORMULATION 
AND OUTPUT FOR BUDGET-CONSTRAINED AND AHP PREFERRED MIX 
MODEL 


SEDs 
I type of combat ships / CG cruisers 
DDG guided missile destroyers 


FFG guided missile frigates 
FSG guided missile corvettes 
PBFA guided missile fast patrol boats / 
P parameters for each ship type 
/C_ cost of ship type I in billions of dollars / 


TABLE 
VALUES (1,P) parameter values for each ship type 
c 


CG 09 
DDG 0.6 

FFG 03 

FSG 02 

PBFA 0.12 ; 


PARAMETER 
NUM (I) _ desired number of ship type I (unattainable) 
OC = 0 
DDG 20 
FFG 20 
FSG 20 
PBFA 20 /; 


SCALARS 
TOTBUDGET /2.6 /; 


INTEGER VARIABLES 


X (I) number of ships I to purchase 
DEVNEG (I) negative deviation from desired number of ships of type I 
DEVPOS (I) positive deviation from desired number of ships of type I 


FREE VARIABLE 
DEVIATION deviation from objective function ; 


EQUATIONS 
BUDGET _ stay within budget 
NSHIPS(I) desired number of ship type I (unattainable) 


OBJECFUN objective function; 


BUDGET.. SUM(I, X(I) * VALUES(I'C’)) =L=TOTBUDGET ; 
qd 


NSHIPS(1).. X(1) - DEVPOS(I) + DEVNEG(I) =E= NUM(I) ; 

OBJECFUN.. 1.9*(DEVNEG( 'CG' ) **2)+ 1.9*DEVNEG (‘CG’) + 
1.4*(DEVNEG( 'DDG' ) **2)+ 1.4*DEVNEG (‘DDG')+ 
1.1*(DEVNEG( 'FFG' ) **2)+ 1.1*DEVNEG (‘FFG')+ 
0.8*(DEVNEG( FSG! ) **2)+ 0.8*DEVNEG (‘FSG')+ 
0.5*(DEVNEG( 'PBFA' ) **2)+0.5*DEVNEG (‘PBFA’) 
=E= DEVIATION ; 


MODEL GP /ALL/: 

SOLVE GP USING RMINLP MINIMIZING DEVIATION: 

B. GAMS QUADRATIC MIXED INTEGER PROBLEM OUTPUT 
SOLVE SUMMARY 


MODEL GP OBJECTIVE DEVIATION 
TYPE? | VIENIEP DIRECTION MINIMIZE 
SOLVER MINOSS5 FROM LINE 77 


**** SOLVER STATUS 1 NORMAL COMPLETION 
*¥*** MODEL STATUS 2 LOCALLY OPTIMAL 
**** OBJECTIVE VALUE 2049.9530 
RESOURCE USAGE, LIMIT 0.110 50000.000 
ITERATION COUNT, LIMIT 4 20000 
EVALUATION ERRORS 0 0 

LOWER LEVEL UPPER MARGINAL 
---- EQU BUDGET -INF 2.600 2.600 -131.469 

LOWER LEVEL UPPER MARGINAL 

---- EQU OBJECFUN , -1.000 
---- VAR X number of ships I to purchase 


LOWER LEVEL UPPER MARGINAL 


CG . 100.000 40.422 
DDG . 100.000 21.482 
FFG : 4.000 100.000 
FSG 4.000 100.000 EES 


PBFA ; 5.000 100.000 EPS 
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---- VAR DEVNEG negative deviation from desired number of ships of type I 


LOWER LEVEL UPPER MARGINAL 


CG 20.000 100.000 
DDG 20.000 100.000 
FFG 16.000 100.000 
HG 16.000 100.000 
PBFA 15.000 100.000 


---- VAR DEVPOS positive deviation from desired number of ships of type I 


LOWER LEVEL UPPER MARGINAL 


CG 100.000 77.900 
DDG 100.000 57.400 
FFG ; 100.000 39.44) 
FSG . : 100.000 26.294 
PBFA 100.000 1/76 


LOWER LEVEL UPPER MARGINAL 


---- VAR DEVIATION -INF 2049.953 +INF 


1 
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